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SUMMARY 


This  report  presents  the  results  of  an  experimental  investiga¬ 
tion  of  rotor  blade  dynamic  stall.  Forces  and  moments  in  two- 
dimensional  flow  on  two  typical  helicopter  rotor  airfoils 
(NACA  0012  and  Vertol  23010-1.58),  in  oscillatory  pitch  and 
vertical  translation,  were  determined  by  measuring  differential 
pressures.  The  Mach  number,  Reynolds  number,  and  ang.’.e  of 
attack  prevailing  in  the  retreating-blade  stall  region  were 
investigated  up  to  typical  first-bending  and  first-torsion  mode 
natural  frequencies. 

Pitch  oscillation  was  found  to  increase  the  angle  of  attack  at 
which  stall  occurred,  with  a  consequent  increase  in  maximum 
normal  force.  Hysteresis  effects  on  the  pitching  moment  re¬ 
sulted  in  negative  aerodynamic  damping  of  oscillations  about 
mean  angles  of  attack  near  the  steady-flow  stall  value  at 
Mach  0.2  and  0.4.  At  Mach  0.6,  little  or  no  such  hysteresis 
appeared.  The  increase  in  usable  angle  of  attack  and  maximum 
normal  force  due  to  airfoil  camber  was  found  to  prevail  under 
nearly  all  oscillatory  conditions  tested,  as  well  as  in  steady 
flow. 

Oscillatory  translation  at  angles  of  attack  near  stall  also 
showed  increased  maximum  normal  force  for  both  airfoils. 
Furthermore,  lag  effects  in  the  stall-unstall  process  were 
noted  which  could  result  in  self-sustaining  blade  bending 
oscillations. 

Additional  tests  were  conducted  with  a  torsionally  flexible 
airfoil  dynamic  system,  tuned  to  frequencies  corresponding  to 
typical  rotor  blade  first-torsion  mode  values,  and  subjected 
to  angle-of-attack  oscillation  at  a  frequency  corresponding  to 
a  helicopter's  1/rev  angle-of-attack  change.  Torsional  oscil¬ 
lations  due  to  stall  were  noted  under  several  conditions. 
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FOREWORD 


The  results  from  the  oscillating  airfoil  tests  are  summarized 
in  this  report.  The  project  was  performed  under  Contract 
DA  44-177-AMC-438{T)  under  the  technical  cognizance  of  Clifton 
G.  Wrestler,  Jr.,  of  the  Aeromechanics  Division  of  USAAVLABS. 

«  The  report  consists  of  two  volumes: 

Volume  I,  Summary  and  Evaluation  of  Results 

Volume  II,  Data  Report 

The  tests  were  conducted  at  the  Commercial  Airplane  Division  of 
The  Boeing  Company  in  the  supersonic  wind  tunnel.  The  assis¬ 
tance  and  cooperation  of  the  Model  Design,  Instrumentation,  and 
Supersonic  Wind  Tunnel  Testing  Groups  are  gratefully  acknowl¬ 
edged. 

Richard  R.  Pruyn  of  the  Vertol  Division  made  significant  con¬ 
tributions  to  this  project  through  his  suggestions  on  model 
design,  instrumentation,  and  data  system  methods. 
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INTRODUCTION 


Contemporary  helicopters  frequently  fly  with  the  rotor  blade 
partially  stalled  in  some  regions  of  the  azimuth.  This  is 
caused  by  the  combined  effects  of  reduced  dynamic  pressure  on 
the  retreating  side  of  the  disc  and  the  requirement  for  roll 
equilibrium.  Furthermore,  in  the  reverse-flow  region,  the 
blade  actually  moves  through  the  air  with  the  trailing  edge 
foremost.  Consequently,  an  increase  in  aircraft  forward  speed 
without  a  corresponding  increase  in  rotor  rpm  implies  en¬ 
largement  of  the  region  where  stall  may  be  a  problem.  In  fact, 
retreating  blade  stall  is  generally  the  factor  which  limits 
helicopter  speed,  regardless  of  the  engine  power  available. 

The  capability  to  predict  the  effects  of  blade  stall  accurately 
is  therefore  a  vital  factor  in  the  successful  development  of 
the  full  speed  potential  of  a  helicopter. 

Stall  prediction  is  much  more  complex  on  helicopter  rotors 
than  on  conventional  aircraft  wings,  because  of  the  hiqhly 
variable  aerodynamic  environment  of  the  rotor  blade  as  it 
moves  around  the  azimuth.  This  environment  is  complex  because: 

1.  The  angle  of  attack  is  influenced  by  blade  motion, 
cyclic  pxteh,  and  the  complex  induced-flow  field 
associated  with  the  blade  vortex  system. 

2.  The  blade  is  practically  never  subjected  to  a  simple 
flow  approaching  j t  at  right  angles  to  its  edges. 
Normally,  there  is  a  substantial  component  of  wind 
velocity  along  the  blade.  This  is  analogous  to  the 
effect  of  sweepback  on  conventional  wings,  which  is 
known  to  have  considerable  influence  on  stall  behavior. 

3.  There  is  the  effect  of  unsteady  flow  on  the  behavior 
of  the  blade  aerodynamic  forces. 

Attempts  to  predict  the  effects  of  blade  stall  on  performance 
by  using  steady-state  airfoil  characteristics  with  stall  and 
Mach  number  accountability  in  a  strip  theory  analysis  have  been 
unsuccessful;  the  onset  of  stall  is  predicted  too  early. 

This  is  shown  in  Figure  1  for  the  H-21  rotor,  which  was  tested 
in  the  40-  by  80-foot  wind  tunnel  at  Ames  Research  Center!. 

The  theory  predicts  stall  at  a  much  lower  value  of  lift  than 
the  wind  tunnel  test  data.  Simple  modifications  to  the  perfor¬ 
mance  theory,  consisting  of  the  deletion  of  stall  in  the  Cl 
data  up  to  an  angle  of  attack  of  20  degrees  (the  drag  and 
pitching  moment  data  are  left  unchanged) ,  have  provided  the 
means  for  predicting  rotor  performance  adequately;  this  is 
shown  in  Figure  1  by  the  dotted  line.  However,  the  angle  of 
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ROTOR  LIFT  FORCE  -  POUNDS  X  1000 


0  1  2  3  4  5 

ROTOR  DRAG  FORCE  -  POUNDS  X  1000 


NOTES:  1.  TEST  CONDITIONS: 

A.  0  75  =  15° 

B.  \t '  =  0.34 


C’  M(  1)  (90)  "  0<  716 


2.  REFERENCE:  NACA  TN  4  367 


3.  THEORY  COMPENSATED  BY  AN 
INCREMENTAL  DRAG  FORCE  DUE 
TO  THE  UNCERTAIN  LEVEL  OF 
WALL  INTERFERENCE  AND  UNKNOWN 
CONTRIBUTION  OF  ROTOR  HUB  DRAG. 


Figure  1.  Rotor  Drag  Correlation  for  H-21  Rotor. 
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attack  for  stall  on  the  rotor  blade  is  somewhere  between 
static  stall  and  the  value  of  20  degrees  which  was  chosen  for 
the  analysis,  and  the  rotor  stall  limit  must  therefore  be 
determined  separately  by  testing  a  similar  model  or  a  full- 
scale  rotor  in  a  wind  tunnel. 

2 

It  has  been  suggested  by  Harris  that  the  increase  in  the 
maximum  C above  the  static  stall  values,  can  be  caused  by 
radial  flow  or  unsteady  aerodynamic  effects.  Radial  flow 
effects  were  included  in  a  performance  analysis,  but  they  did 
not  increase  Cl  sufficiently  to  correlate  theory  with  test  data; 
stall  was  still  predicted  too  early. 

3 

Half man's  test  data  on  a  pitching  airfoil  show  that  the  maxi¬ 
mum  lift  experienced  by  an  airfoil  oscillating  through  stall 
can  be  substantially  higher  than  the  two-dimensional  static 
value.  Dynamic  effects,  then,  are  the  most  likely  reason  for 
the  stall  delay  exhibited  by  a  helicopter  rotor  blade. 

Another  problem  is  also  important  for  safety  reasons.  Flight 
test  experience  has  indicated  that  large  blade  torsional  oscil¬ 
lations  can  be  encountered  before  severe  vibration  warns  the 
pilot  of  stall,  especially  under  high  weight  conditions.  More¬ 
over,  this  can  happen  before  the  speed  limit  due  to  maximum 
available  engine  power  is  reached. 

This  phenomenon  involves  interaction  of  the  blade’s  torsional 
elastic  behavior  with  its  aerodynamic  pitching  moment  (Cm) . 

Under  the  high-a,  unsteady-flow  conditions  in  question,  Cm  be¬ 
haves  in  a  highly  nonlinear  way.  Therefore,  to  provide  a 
basis  for  predicting  the  onset  and  severity  of  such  oscilla¬ 
tions,  experimental  Cm  data  for  those  conditions  are  needed. 

A  literature  search  (see  SELECTED  BIBLIOGRAPHY)  reveals  that 
even  though  many  wind  tunnel  tests  have  been  performed  to 
evaluate  the  effects  of  oscillatory  motion  on  airfoil  loads, 
most  of  them  have  been  aimed  at  establishing  effects  on  air¬ 
plane  wings;  thus,  they  were  carried  out  under  test  conditions 
and  for  model  configurations  which  are  net  directly  applicable 
to  helicopter  rotor  blades. 

The  only  such  data  directly  applicable  (pitching  about  quarter- 
chord,  at  realistic  Mach  and  Reynolds  numbers)  were  obtained 
by  Carta.  His  data,  however,  remain  unavailable  except  for 
some  average  damping  curves  and  the  few  normal  force  and  pitch¬ 
ing  moment  curves  published'*  >  5 .  Halfman's  data,  already  noted, 
did  cover  both  pitching  and  translation  motion  of  a  typical 
helicopter  airfoil.  However,  the  Mach  number  range  was  not 
high  enough,  and  the  pitch  axis  was  too  far  aft. 

Model  and  full-scale  rotor  tests  are  another  source  of  lift 
and  pitching  moment  data.  However,  such  tests  are  highly 
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three-dimensional,  and  it  is  difficult  to  separate  the 
contributions  of  radial  flow  from  'those  of  unsteady  flow. 
Furthermore,  the  angle  of  attack  is  a  very  uncertain  quantity. 
These  data  are  still  valuable  in  that  they  can  provide  an 
ultimate  check  on  the  analytical  methods  which  must  be  developed 
to  use  nonsteady  two-dimensional  airfoil  data. 

To  provide  this  needed  capability  to  predict  stall  under 
dynamic  conditions,  both  for  rotor  loads  and  for  performance 
analyses,  the  program  to  wind  tunnel  test  oscillating  airfoils 
reported  herein  was  undertaken.  Its  objective  was  to  provide 
a  comprehensive  framework  of  experimental  data  around  which 
theories  of  dynamic  stall  can  be  developed. 
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SELECTION  OF  TEST  CONDITIONS  AND  VARIABLES 


The  model  characteristics  and  the  operating  conditions  re¬ 
quired  for  a  comprehensive  series  of  tests  designed  to  isolate 
and  clarify  the  influence  of  unsteady  flow  on  rotor  blade 
stall  are  reviewed  in  this  section. 

The  five  items  to  be  defined  are: 

1.  Mode  of  testing  (i.e.,  two-  versus  three-dimensional 
flow) 

2.  Airfoil  configuration 

3.  Model  scaling 

4 .  Airfoil  motion 
TEST  MODE 


The  purpose  of  this  test  program  was  to  isolate  the  effects  of 
time-varying  angle  of  attack  on  rotor  stall.  A  two- 
dimensional  flow  has  all  the  main  characteristics  which  will 
create  the  pressure  distribution  traits  responsible  for 
boundary  layer  separation.  A  two-dimensional  wind  tunnel  also 
provides  very  close  control  of  airfoil  angle  and  motion  with 
respect  to  the  relative  wind;  therefore,  such  a  facility  was 
a  natural  choice. 

The  absence  of  a  centrifugal  force  field  and  of  a  fluctuating 
spanwise  component  of  free-stream  flow  can  have  only  secondary 
effects  on  the  results  below  and  at  stall.  Their  influence  on 
the  blade  surface  iressures  arises  through  changes  in  boundary 
layer  thickness  due  to  spanwise  variations  in  pressure  or 
velocity.  These  changes  are  small  compared  to  the  blade  sur¬ 
face  motions  of  interest  here. 

AIRFOIL  CONFIGURATION 


To  establish  a  link  to  previous  work  and  to  provide  a  data  base 
against  which  more  advanced  airfoil  sections  can  be  measured, 
the  NACA  0012  (modified)  was  the  first  airfoil  chosen  for 
testing.  The  modifications  consist  of  a  symmetrical  Jeading- 
edge  fairing  (corresponding  to  a  typical  anti-icing  boot)  and 
a  flat  sheet  metal  trailing-edge  extension.  The  resulting 
configuration  is  identical  to  the  rotor  blade  section  used  on 
the  CH-47A  helicopter.  An  advanced  airfoil,  of  the  type  used 
on  the  CH-47B,  was  chosen  for  comparison.  This  airfoil  has  a 
camber  line  similar  to  the  NACA  230  and  a  thickness  ratio  of 
about  10  percent.  A  cuspate  trailing  edge  similar  to  that  of 
the  first  airfoil  is  used.  The  second  airfoil  was  also  tested 
with  a  3-degree  upward  deflection  of  the  trailing  edge,  which 
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Airfoil 


flight  test  experience  has  shown  to  be  desirable, 
contours  and  coordinates  are  shown  in  Figure  2. 

MODEL  SCALING 


Three  basic  parameters  govern  the  scaling  of  the  model.  These 
are  the  Mach  number  (M)  ,  Reynolds  number  (R)  ,  and  reduced 
frequency  (k) . 
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For  full  aerodynamic  similarity,  the  Mach  number  and  the 
Reynolds  number  for  the  model  and  the  full-scale  rotor  blade 
must  be  identical. 


(R) MODEL 

(VC) MODEL 

(v)Fl!LL  SCALE 

(r)full  SCALE 

(VC}FULL  SCALE 

(v) MODEL 

(M) MODEL 

(V)MODEL 

(a)FULL  SCALE 

TmJ 

'  'FULL  SCALE 

(V) 

v  'FULL  SCALE 

(a) MODEL 

Equations  (4)  and  (5)  are  satisfied  by  testing  at  full-scale 
velocity  in  a  pressurized  wind  tunnel  with  air  as  the  testing 
n  dium.  Since  the  tunnel  size  dictated  a  model  chord  of  27 
percent  of  full-scale  size,  the  test  section  total  pressure 
and  the  kinematic  viscosity  were  adjusted  to  provide  full- 
scale  Reynolds  numbers. 

Figure  3  shows  the  Mach  and  Reynolds  numbers  used  for  these 
tests.  Figure  4  shows  the  advancing  and  retreating  tip  Mach 
number  conditions  typical  of  current  and  advanced-design  heli¬ 
copter  rotors.  Mach  numbers  from  0.3  to  0.5  are  of  special 
interest  at  the  tip  of  the  retreating  blade  near  the  270-degree 
azimuth.  However,  conditions  at  other  points  (highest  angles 
of  attack  are  often  reached  at  approximately  300  degrees)  re¬ 
quire  the  somewhat  wider  band  of  0.2  to  0.6  Mach  number  chosen 
for  these  tests. 

For  dynamic  similarity,  the  reduced  frequency  for  the  model  and 
full-scale  blade  must  be  identical.  Since  the  velocity  V  was 
determined  by  aerodynamic  considerations,  the  requirement  for 
dynamic  similarity  can  be  satisfied  as  follows: 
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Figure  2.  Test  Airfoil  Coordinates. 
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Figure  3.  Reynolds  Number  Variation  With  Mach  Number 
for  Wind  Tunnel. 
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Figure  4.  Present  and  Envisioned  Helicopter  Operating  Region- 


(k)MODEL  _  (fc)MODEL 

( f  c )  < 

K  FULL  SCALE  '  ;  FULL  SCALE 

This  was  achieved  by  increasing  the  model  oscillation  frequency 
in  proportion  to  the  decrease  in  model  chord.  The  testing 
frequencies  and  amplitudes  of  motion  are  discussed  in  the  next 
section . 

AIRFOIL  MOTIONS 

Before  considering  what  kinematic  characteristics  will  be  re¬ 
quired  of  the  model,  the  following  should  be  noted  concerning 
the  question  of  frequency  scaling  between  the  airfoil  model 
and  the  actual  rotor. 

Aerodynamic  similarity  requires  that  the  time  to  execute  a 
given  motion,  expressed  as  a  ratio  to  the  time  for  the  air- 
stream  to  move  one  semichord,  should  be  the  same  for  the  model 
as  for  the  rotor  blade.  It  can  easily  be  shown  that  if  this 
quantity  is  matched  (along  with  the  Mach  and  Reynolds  numbers)  , 
then  complete  two-dimensional  aerodynamic  similarity  is 
obtained.  For  periodic  motions,  this  characteristic  time  is 
used  to  nondimensionalize  the  frequency  to  obtain  k,  the 
reduced  frequency  (sometimes  called  the  Strouhal  number)  . 

To  investigate  the  required  airfoil  kinematics,  consider  first 
the  rotor  blade  relative-wind  vector  diagram  sketched  below. 


The  following  time  variations  appear  in  the  parameters  shown: 

1.  uT 


a.  Varies  sinusoidally  at  1/rev  because  of  aircraft 
forward  speed 


b.  Varies  at  higher  frequencies  (but  very  low 
amplitude)  because  of  chord  bending 

c.  Varies  irregularly  because  of  vortex  effects 


a.  Has  no  1/rev  variation  since  axes  are  referred  to 
plane  of  no  flapping 

b.  Varies  at  frequencies  of  2/rev  and  higher  because 
of  flap  bending  at  appreciable  amplitude 

c.  Varies  irregularly  because  of  vortex  effects 
3’  A 

a.  Varies  at  1/rev  because  of  cyclic  pitch  and/or 
flapping  motion  (accepting  the  classical  equiva¬ 
lence  between  flapping  and  feathering) 

b.  Varies  at  higher  frequencies  (4/rev  or  higher) 
because  of  blade  torsional  motion 

It  should  be  noted  in  regard  to  Oa  that  typical  hub  and  rotor 
blade  designs  have  both  the  pitch  bearing  axis  and  the  shear 
center  near,  or  slightly  forward  of,  the  quarter-chord  point; 
therefore,  that  location  on  the  airfoil  chord  line  was  chosen 
as  the  center  of  rotation  for  pitch  oscillations. 

For  the  rotor  blade,  U,  <f> ,  and  a  are  dependent  variables  which 
can  be  deduced  from  the  three  quantities  already  discussed. 

For  the  wind  tunnel  model, the  independent  variables  can  also 
be  taken  as  U-p,  Up,  and  ^a*  (In  the  discussion  below,  the 
corresponding  wind  tunnt  variables  will  be  denoted  by  IT  ,  Up , 
and  oA  to  avoid  confusion.)  U-p  Up,  and  Oa  are  subject  to 
certain  restrictions  because  of  mechanical  or  operational 
limitations  of  the  model  or  test  facility;  i.e., 

1.  UT  may  not  oscillate. 

2.  Up  is  limited  to  sinusoidal  variations  of  amplitude 
fAh  about  a  mean  value  of  zero.  Furthermore,  as  f 
increases,  the  available  Ah  will  be  limited  by 
strength  and  vibration  considerations,  and  Ah  cannot 
exceed  0.31  chord  in  any  case. 

3.  Oa  may  oscillate  over  a  large  range  of  amplitude  and 
frequency  at  any  desired  mean  value;  however,  0a  and 
Up  cannot  oscillate  at  the  same  time. 
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Within  these  constraints,  which  of  the  time-varying  parameters 
listed  above  can  be  studied  meaningfully?  Consider  them  in 
the  order  given: 

1.  Direct  simulation  of  Ut  variation  by  Ut  is  obviously 
ruled  out.  However,  it  is  conceivable  that  a  vari¬ 
ations  due  to  changes  in  Ut  (item  la)  at  fixed  Up 
might  be  simulated  by  0^  oscillations  at  constant  Ut- 
(A  substantial  average  negative  Up  must  prevail  on 

a  rotor  producing  propulsive  force  at  high  speed.) 

This  possibility  was  explored  analytically  using 
unsteady-flow  potential  theory.  It  was  concludrd  that 
such  a  simulation  would  not  be  meaningful,  because  of 
time  lag  effects  due  to  the  shedding  of  counter¬ 
vortexes  . 

2.  Up  variation  due  to  flap  bending  (item  2b)  can  be 
simulated  and  is  of  interest.  Appropriate  frequencies 
and  amplitudes  are  discussed  in  a  following  paragraph 
(Blade  Elastic  Motion  Study) .  Item  2c  was  considered 
to  be  beyond  the  scope  of  the  work  undertaken  in  this 
effort . 

3.  The  1/rev  variation  of  0A_(item  3a)  can  easily  be 
simulated  by  appropriate  9^  oscillations.  Dynamic 
lift  curve  behavior  above  static  stall,  under  these 
conditions,  was  one  of  the  two  main  concerns  prompting 
this  investigation.  Item  3b,  higher  frequency  0 
variations,  can  also  be  studied  by  0A  motions;  this  is 
the  second  principal  reason  for  undertaking  this  work. 
Frequencies,  amplitudes,  and  other  aspects  of  the 
motion  are  also  discussed  in  the  next  major  paragraph. 

It  will  not  be  necessary  again  to  discuss  relative  wind,  angle 
of  attack,  or  other  elements  of  the  airfoil  model  as  distinct 
from  those  of  the  rotor  blade  in  the  same  context.  Therefore, 
the  notation  used  in  the  remainder  of  this  report  will  conform 
to  the  usual  wind  tunnel  practice/  i.e.,  a  (instead  of  0^) 
will  refer  to  the  angle  between  airfoil  chord  line  and  the 
wind  tunnel  centerline,  and  V  (instead  of  UT)  will  refer  to 
test  section  wind  speed. 

BLADE  ELASTIC  MOTION  STUDY 


A  computer  analysis  was  made  of  rotor  blade  aeroelastic 
motions  for  a  variety  of  blade  configurations  and  flight  con¬ 
ditions.  The  ranges  of  amplitude  and  frequency  of  those 
results  were  used  as  guides  in  defining  test  conditions. 
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The  results  were  as  follows: 

1.  Pitching  motion  about  the  airfoil  quarter-chord  at 
amplitudes  of  2.5  degrees  and  5  degrees  and  fre¬ 
quencies  corresponding  to  the  range  from  4/rev  to 
6/rev  were  found  to  be  sufficient  to  cover  the  effects 
of  blade  torsional  oscillations.  To  permit  clear 
definition  of  trends  with  frequency,  it  was  also  found 
necessary  to  conduct  tests  of  pitching  motions  at  2/rev 
and  3/rev.  The  higher  frequency  pitching  tests, 
together  with  1/rev  oscillations  at  amplitudes  up  to 
7.5  degrees,  provided  thorough  coverage  of  the  ranges 

of  interest  in  rotor  technology. 

2.  Plunging  motions  of  0.3  and  0.45  semichord  amplitude, 
at  a  frequency  corresponding  to  2/rev,  were  found  to  be 
sufficient  to  cover  the  effects  of  blade  bending.  To 
establish  trends  with  frequency,  tests  were  also  found 
to  be  desirable  at  1/rev. 

STALL  OF  A  SIMULATED  TORSIQNALLY  FLEXIBLE  BLADE 

When  stall  occurs,  the  actual  pitching  motion  of  a  rotor  blade 
element  will  generally  be  much  more  complex  than  a  simple 
harmonic  oscillation.  Becau.se  rotor  blades  are  usually  very 
limber  in  torsion,  the  large  nose-down  pitching  moment  normally 
associated  with  stall  can  cause  a  rapid  and  substantial  nega¬ 
tive  elastic  twist.  This  leads  to  loss  of  lift  due  to  reduced 
angle  of  attack  and  other  complications.  (Harris  and  Pruyn 
discuss  this  in  their  paper^.)  Under  some  circumstances,  self- 
sustaining  limit-cycle  oscillations  (or  stall  flutter)  may 
occur,  as  discussed  by  Carta'*. 

Motions  of  two  distinct  types  appear  to  be  involved:  (1)  the 
1/rev  sinusoidal  variation  due  to  cyclic  pitch  control,  and 
(2)  torsional  oscillations  which  if  present  would  be  expected 
to  occur  at  a  frequency  close  to  the  first  torsional  natural 
frequency  of  the  blade  (in  the  range  from  4/rev  to  6/rev). 

These  motions  are  superposed  on  the  1/rev  motion.  Their 
amplitude  and  their  persistence  with  time  will  depend  on  the 
blade  dynamic  characteristics  and  on  the  amplitude  and  mean 
value  of  the  1/rev  background  motion. 

In  order  to  show  the  combined  effects  of  the  two  kinds  of 
motion,  it  was  decided  to  conduct  tests  of  a  tuned  system,  in 
which  the  airfoil  is  mounted  on  a  torsion  spring  through  which 
the  background  1/rev  motion  is  transmitted.  The  spring  stiff¬ 
ness  was  varied  to  give  a  range  of  natural  frequencies  from 
four  to  six  times  the  background  frequency. 
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The  data  from  these  tests  also  serve  a  second  purpose: 
represent  a  comparatively  simple  dynamic  system  on  whic 
mathematical  models  constructed  from  the  data  obtained 
forced  sinusoidal  motion  may  be  tried  out. 
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TEST  FACILITIES,  DATA  RECORDING, 
AND  REDUCTION  SYSTEMS 


A  detailed  description  of  the  test  facility,  apparatus,  and  data 
reduction  system  is  presented  in  Volume  II  of  this  report.  A 
short  summary  is  presented  here. 

TEST  FACILITIES 

The  variable  density,  two-dimensional  l-by-3-foot  test  section 
of  the  Boeing  4-by-4-foot  supersonic  wind  tunnel  was  used  for 
these  tests.  The  tunnel  total  pressure  and  the  resulting  air 
density  were  varied  to  simulate  full-  and  half-scale  Reynolds 
numbers  for  a  typical  transport  helicopter  rotor  blade. 

Pitch  Mechanism 


The  pitch  oscillating  mechanism  is  shown  in  Figure  5.  A  fly¬ 
wheel  with  an  eccentrical ly  mounted,  interchangeable  cam 
drives  a  crank  which  is  connected  to  the  airfoil  by  a  torsion 
spring  for  the  tuned  system,  or  by  an  adapter  for  direct 
drive.  A  hydraulic  motor  supplies  the  power.  An  eddy-current 
damper  was  used  with  the  tuned  system  for  safety. 

Vertical  Translation 


A  system  similar  to  the  pitch  mechanism  was  used  for  the  trans¬ 
lation  tests.  Airfoil  endplat.es  were  required  with  the  trans¬ 
lating  mechanism  to  cover  the  slots  in  the  tunnel  walls. 

These  end  plates  caused  flow  separation  at  the  tunnel  wall  which 
was  eliminated  by  using  wall  boundary  layer  control. 

Models  and  Instrumentation 


The  model  airfoils  were  constructed  of  a  steel  center  spar  with 
sheet  metal  leading  and  trailing  edges.  The  span  of  each  wing 
was  nominally  12  inches  and  the  chord  was  6.38  inches.  Four¬ 
teen  miniature  differential-pressure  transducers  were  used  to 
obtain  the  chordwise  pressure  distribution.  The  natural 
frequency  of  the  installed  transducers  was  at  least  1030  Hertz 
(Hz)  with  a  damping  factor  of  0.45. 

DATA  RECORDING  SYSTEM 

All  data,  both  airfoil  pressure  and  tunnel  test  condition 
information,  were  recorded  on  two  FM  wideband  tape  recorders. 

A  time  code  and  a  l/'rev  pulse  recorded  on  both  tapes  provided 
tape  synchronization.  In  addition,  all  parameters  were  record¬ 
ed  on  oscillograph  tape  for  instant  data  monitoring.  Included 
were  Cn  and  CM,  obtained  by  integrating  the  pressures  with  an 
analog  computer.  The  data  recording  flow  path  is  shown  in 
Figure  6. 


15 


flywheel 


CRANK 


TURNTABLE  (SETS  a.) 


PRESSURE  TRANSDUCER  PORTS 


HYDRAULIC 


ADJUSTABLE  ECCENTRIC 


DAMPER  (SCHEMATIC) 


ADAPTER  FOR  DIRECT  DRIVE 


Figure  5.  Pitch  Oscillating  Mechanism, 
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Figure  6.  Data  Recording  System. 
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DATA  REDUCTION  SYSTEM 

Nearly  800  individual  test  points  were  taken  in  a  series  of 
four  tests.  The  total  number  of  individual  data  readings 
was  close  to  10  million,  since  the  data  for  10  cycles  were 
averaged  and  approximately  40  points  were  read  per  cycle.  To 
handle  this  large  volume  of  data,  the  recorded  analog  informa¬ 
tion  was  digitized  and  processed  on  a  digital  computer,  using 
the  basic  flight  test  data  reduction  programs  which  we;re  used 
by  Pruyn^  for  the  "Tandem  Rotor  Airloads  Measurement  Program". 
The  pressure  data  were  integrated  to  obtain  time  histories  of 
CN,  CM,  and  W,  the  work  or  damping  parameter,  and  the  average 
cycle  damping .  The  flight  test  data  plotting  system,  which  is 
integral  with  the  data  reduction  system,  was  used  to  plot 
approximately  80  percent  of  the  CN  and  CM  versus  a  data. 
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EXPERIMENTAL  RESULTS 


FORCED  PITCHING  OSCILLATION 

The  forced  pitching  oscillation  data  points  cover  three  airfoil 
configurations  (symmetrical,  cambered,  and  cambered  with  reflex) 
for  mean  angles  of  attack  from  0  to  25  degrees,  amplitudes  from 
2.5  to  7.5  degrees,  Mach  numbers  from  0.2  to  0.6,  and  reduced 
frequencies  from  0.042  to  0.720. 

The  effects  of  the  flow  and  oscillation  parameters  on  the  CN 
and  CM  versus  a  behavior  of  the  Vertol  23010-1.58  airfoil  are 
presented  in  Figures  7  through  12.  (The  general  features  and 
trends  are  very  similar  for  both  the  cambered  and  the  symmet¬ 
rical  sections.)  The  significant  differences  in  behavior  of  the 
airfoils  are  then  discussed,  and  data  for  the  NACA  0012  profile 
are  presented  which  illustrate  those  differences. 

Effect  of  Mean  Angle  of  Attack 

Figure  7  presents  a  set  of  C^  and  Cm  versus  a  traces  for  5- 
degree  amplitude  oscillations  of  the  23010-1.58  airfoil  at 
M  =  0.4  and  k  =  0.06.  A  sequence  of  increasing  o0's  is  shown 
and  compared  with  steady-a  data. 

The  first  pair  of  traces,  a0  =  7.33  degrees,  shows  the  charac¬ 
teristic  elliptical  shape  centered  on  the  steady-a  line  of 
stall-free  flow.  The  displacements  from  steady-a  data  are  con¬ 
sistent  with  potential  theory.  The  next  pair,  at  ac  =  9.63 
degrees,  shows  the  first  signs  of  stall  in  the  figure-eight 
shapes  of  both  Cm  and  Cm-  At  the  high-a  end  of  the  cycle,  loops 
are  formed  in  a  clockwise  sense. 

For  Cm,  the  area  enclosed  by  the  trace  and  the  sense  of  motion 
around  the  loops  have  important  physical  significance.  The  net 
work  done  by  the  airfoil  on  the  surrounding  air  is  proportional 
to  the  integral 

W  =  i  Cu  da.  (7) 

This  integral  is  proportional  to  the  area  enclosed  by  the  trace 
and  is  positive  for  a  counterclockwise  (CCW)  circuit.  If  the 
circuit  encloses  a  substantial  area  in  a  clockwise  (CW)  sense, 
the  contribution  of  that  area  is  negative;  i.e.,  it  represents 
energy  extracted  from  the  airstream  by  the  airfoil.  Net  energy 
extraction  in  a  cycle  (negative  damping)  implies  that  the  rotor 
blade  oscillation  in  which  it  occurred  would  tend  to  increase 
in  amplitude.  This  is  precisely  the  condition  for  flutter. 

The  third  pair,  at  a0  =  12.22  degrees,  shows  substantial  stall 
effects;  e.g.,  large  departures  from  zero  Cm  and  considerable 
loss  of  lift.  Some  aspects  of  these  traces  worthy  of  special 
note  are: 
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1.  The  angles  of  attack  at  which  CM  and  CN  break  are 
clearly  higher  than  the  steady-a  values.  CN  breaks 
at  16  degrees  and  CM  at  lb  degrees,  where  both  break 
at  14  degrees  in  steady  flow. 

2.  A  large  CW  loop  is  evident  in  the  middle  of  the 
trace.  The  net  damping  for  this  cycle  is  negative. 
However,  a  CCW  loop  has  appeared  at  the  high-a  end 
of  the  cycle,  indicating  that  as  aQ  increases,  there 
will  be  a  return  to  more  stable  conditions. 

3.  There  is  a  substantial  lag  in  tne  reestablishment  of 
attached  flow  after  a  has  dropped  below  the  steady-state 
flow  stall  level. 

The  traces  for  a0  =  14.92  degrees  show  further  development  of 
the  same  features.  The  CCW  loop  at  high  a  has  now  become  large 
enough  to  balance  the  big  CW  loop  in  the  center,  and  the  damp¬ 
ing  is  only  slightly  negative.  For  a0  =  17.51  degrees,  the  net 
damping  is  again  positive.  It  is  worth  noting  that  the  angles 
for  the  Cn  and  Cm  break  have  remained  at  17  and  15  degrees, 
respectively,  for  the  last  three  conditions  discussed. 

The  final  set,  at  a0  =  24.57  oegrees,  shows  a  condition  of 
fully  developed  stall  at  all  times.  The  mechanism  of  flow 
separation  and  reattachment  no  longer  creates  the  loop  shape 
required  for  feeding  energy  into  the  airfoil  dynamic  system. 

Effect  of  Pitchirc_  ri^quencv 

Figure  8  shows  C^  and  Cm  traces  for  the  23010-1.58  airfoil  at 
a  sequence  of  increasing  frequencies  for  M  =  0.4,  a0  =  12.5 
degrees,  and  ha  =  5  degrees.  (The  angles  are  nominal  values; 
ha  actually  increases  from  4.83  degrees  to  5.65  degrees  as  k 
goes  from  0.062  to  0.355,  because  of  the  dynamics  of  the 
oscillating  mechanism) .  The  a0  chosen  here  is  of  particular 
significance  because  it  is  in  the  region  of  highest  negative 
damping. 

The  most  significant  trend  with  frequency  is  the  collapse  of 
the  Cn  hysteresis  loop  into  a  thin  ellipse.  The  resulting 
trace  runs  parallel  to  and  slightly  below  the  static  C^g  versus 
a  curve  because  Cn  is  lower  during  the  decreasing-a  part  of 
the  cycle  than  it  is  during  the  increasing-a  part,  contrary 
to  theory. 

The  flattening  of  the  loop  may  be  explained  as  the  result  of 
two  different  phenomena.  The  first  is  the  postponement  of  the 
onset  of  stall  to  the  point  where  a  stops  increasing,  evident 
at  k  =  0.124  and  above.  This  is  reasonable  in  view  of  the 
reduction  of  the  apparent  angle  of  attack  at  the  leading  edge 
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because  of  the  upward  motion  of  the  airfoil.  (Effectively, 
pitching  motion  applies  camber  to  the  boundary  conditions 
governing  the  flow  about  the  airfoil,  in  direct  proportion  to 
pitch  rate.)  The  second  phenomenon  is  an  apparent  upper  limit 
to  the  time  rate-of-change  of  lift. 

The  CN  curve  is  displaced  to  a  lower  level  because  stall  re¬ 
covery  is  not  complete  until  the  second  half  of  the  increasing- 
■  i  cycle  is  reached.  A  theoretical  trace  for  stall-free  flow 
is  shown  on  the  plot  (Figure  8)  for  k  =  0.355  to  confirm  the 
above  point . 

The  CM ' s  also  show  increasing  postponement  of  stall  as  frequency 
goes  up.  The  Cm  break  always  occurs  earlier  in  cycle  time  than 
the  Cfj  break.  This  tendency  of  moment  stall  to  respond  less  to 
pitch  rate  than  does  lift  stall  is  apparently  universal. 

As  the  a  for  moment  stall  goes  up,  the  dynamic  Cm  overshoots  the 
static-stall  Cm  level  by  an  increasing  margin.  Ultimately, 
however,  an  upper  bound  to  the  time  rate-of-change  again  appears, 
moderating  the  overshoot  and  postponing  the  return  to 
potential-flow  Cm  levels.  Therefore,  the  CW  loop  in  the  Cm 
trace  grows  with  increasing  k  until  all  portions  of  the  cycle 
contribute  to  negative  damping.  At  the  highest  k's,  however, 
the  rate-of-change  limit  discussed  above  begins  to  reduce  the 
loop  size. 

The  potential- f low  Cm  ellipse  has  been  aaded  to  the  k  =  0.355 
CM  graph.  In  this  case,  the  test  results  follow  closely  the 
theoretical  lino  over  the  entire  upward-pitching  half  cycle. 

Since  the  behavior  of  the  chordwise  pressure  distribution  (to 
be  discussed  later)  indicates  that  attached  flow  prevails  only 
over  the  increasing-a  part  of  the  cycle,  the  agreement  with 
theory  of  that  part  of  the  moment  trace  can  be  expected. 

Figure  9  shows  Cjg  and  Cm  as  functions  of  airfoil  pitch  refer¬ 
ence  angle.  (These  curves  may  be  considered  as  functions  of 
time,  but  note  that  the  scale  varies  inversely  with  the 
frequency.)  The  concepts  of  lift  and  moment  stall  discussed 
above  are  better  illustrated  in  this  type  of  plot.  It  is  also 
significant  that,  although  the  C^  versus  0  behavior  approaches 
the  sinusoidal  behavior  of  unstalled  flow,  at  high  frequencies 
the  Cm  still  shows  stall  effects  up  to  the  highest  frequency 
of  this  test. 

Effect  of  Pitch  Amplitude  Maximum  a  Fixed 

The  Cn  and  CM  traces  for  three  different  amplitudes  of  oscil¬ 
lation  at  M  =  0.4  and  k  =  0.124  are  compared  in  Figure  10. 

Since  the  maximum  a  and  the  frequency  of  oscillation  are  the 
same  for  each  case,  the  differences  in  C^  and  CM  behavior 
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above  a  =  15  degrees  are  caused  by  pitch  rate  effects  (da/dt 
being  proportional  to  4a)  .  CjgMAX  increases  significantly  with 
increased  du/dt.  Pitching  moment  break  is  also  delayed  by  the 
higher  da/dt  of  the  large-amplitude  oscillation.  In  addition, 
since  the  2.5-degree  oscillation  does  not  extend  into  the  a 
region  below  static  stall,  nearly  the  whole  moment  curve  con¬ 
tributes  to  negative  damping.  For  the  larger  amplitude  oscil¬ 
lation  both  positive  and  negative  damping  are  present,  result¬ 
ing  in  a  net  positive  damping  for  Aa  of  7.5  degrees  and  a 
negative  damping  for  Aa  of  5  degrees. 

Effect  of  Mach  Number 


Figures  11  and  12  show  sequences  of  CM  and  CN  traces  at 
f  =  16  Hertz  (scale  1/rev)  and  Aa  =  5  degrees  for  increasing 
a0  at  M  =  0.2  and  M  =  0.6.  They  should  be  considered  along 
with  the  corresponding  data  for  M  =  0.4,  already  presented 
in  Figure  7. 

Together,  these  three  sequences  illustrate  airfoil  response  to 
the  1/rev  a  variation  of  a  rotor  blade  at  the  Mach  number 
range  appropriate  to  the  retreating  blade.  The  questions  of 
greatest  concern  in  that  response  are: 

1.  Under  what  conditions  does  Cjg  show  abrupt  and/or  sub¬ 
stantial  reductions  from  its  potential-flow  level? 

2.  Under  what  conditions  does  show  large  deviations 
from  zero,  and  how  big  are  the  deviations? 

At  M  =  0.2,  the  trace  for  a0  =  12.28  degrees  shows  no  evidence 
of  stall,  either  in  Cn  or  in  C^.  At  this  point,  large  stall 
effects  were  evident  at  M  =  0.4.  Even  at  a0  =  14.97  degrees, 
an  a  of  20  degrees  is  reached  before  stall,  which  is  later 
precipitated  by  the  adverse  effect  of  nose-down  pitch  rate 
(effective  inverse  camber).  As  usual,  C^  begins  to  decrease 
before  CN,  but  it  does  not  overshoot  the  static  level,  and  it 
varies  at  a  moderate  rate. 

At  a o  =  17.53  degrees,  more  spectacular  events  occur  when  a 
reaches  its  maximum.  The  drops  very  rapidly  to  -0.32, 
three  times  the  static-stall  value.  Its  recovery  is  fast,  but 
this  impulsive  loading  could  evoke  a  strong  torsional  response 
in  a  rotor  blade.  The  Cn  begins  its  maneuver  by  going  up, 
not  down.  The  reason  for  such  unusual  behavior  is  probably 
connected  with  the  fact  that  potential-flow  conditions  had 
been  carried  to  an  exceptionally  high  a  before  breakdown 
occurred.  The  pressure  distribution  for  this  condition  will 
be  examined  later. 
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The  potential-flow  lift  level  is  not  achieved  at  any  point  in 
the  cycle  for  the  a0  =  19.78°Cn  trace.  This  accounts  for  the 
comparatively  quiescent  behavior  of  Cj^. 

At  M  =  0.6,  a  completely  different  pattern  prevails.  The  whole 
set  of  Cn  and  CM  traces  is  characterized  by  minimal  deviations 
from  the  static  values  over  the  entire  range  of  a0's  tested, 
because  of  the  very  low  reduced  frequency  (17  Hertz  corresponds 
to  k  =  0.04  at  this  speed,  as  opposed  to  0.12  at  M  -  0.2)  . 
Furthermore,  the  static  CN  and  CM  curves  change  gradually  with 
a,  and  they  are  free  from  the  discontinuities  evident  at  lower 
Mach  numbers.  The  reasons  for  the  continuous  changes  through 
stall  at  M  =  0 . 6  for  both  airfoils  will  be  discussed  in  a  later 
section . 

Effects  of  Reynolds  Number 

The  Reynolds  number  was  decreased  to  a  value  corresponding  to 
half  scale  (60  percent  for  M  =  0.2)  for  several  test  conditions 
by  reducing  the  total  pressure  in  the  test  section.  No  curves 
are  shown,  because  no  noticeable  differences  between  half-scale 
and  full-scale  data  were  found.  There  was  one  exception:  at 
M  =  0.2,  differences  were  observed  when  the  tunnel  was  operating 
at  its  low-speed  limit  and  the  angle  of  attack  was  doubtful; 
however,  these  differences  are  not  considered  to  be  meaningful. 

Effect  of  Airfoil  Symmetry 

Figure  13  presents  a  sequence  of  traces  for  the  NACA  0012 
(modified)  symmetrical  airfoil  oscillating  +5  degrees  about 
a  o  =  12.3  degrees.  The  Mach  number  is  0.4  and  the  frequency 

varies  from  k  =  0.065  to  k  =  0.371.  This  sequence  can  be  com¬ 
pared  directly  to  that  given  in  Figure  8  for  the  cambered 
airfoil . 

At  the  lowest  frequency  (equivalent  to  1/rev)  ,  the  traces  show 
substantially  earlier  stall  than  was  evident  for  the  23010-1.58 
airfoil.  In  fact,  they  strongly  resemble  the  a0  =  15-degree 
curves  for  that  section  (Figure  8,  fourth  pair). 

As  the  frequency  goes  up,  the  loop  tends  to  flatten,  much  as 
it  did  for  the  23010-1.58  airfoil.  At  the  same  time,  however, 
significantly  larger  0M  excursions  are  evident.  A  minimum  CM 
of  -0.3,  which  does  not  decrease  with  further  increases  in 
frequency,  is  reached  near  k  =  0.2.  Moreover,  recovery  from 
this  large  negative  value  occurs  later  in  the  cycle  as  the 
frequency  is  increased.  Thus,  the  CW  area  in  the  loop,  and 
therefore,  the  negative  damping,  is  greater  for  the  (Symmetrical 
than  for  the  cambered  airfoil. 
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Figure  14  shows  the  influence  of  Mach  number  on  the  stalling 
response  of  the  symmetrical  airfoil  at  f  =  16  Hertz  (equivalent 
to  1/rev).  At  M  =  0.2  and  a0  =  12.3  degrees,  large  stall  effects 
are  evident  in  both  CN  and  CM ,  whereas  the  cambered  section 
showed  essentially  perfect  potential-flow  behavior  under  the 
same  conditions.  (See  Figure  11,  second  pair  of  traces).  At 
M  =  0.6,  the  dynamic  effects  have  essentially  disappeared,  and 
both  Cn  and  Cm  remain  close  to  the  static  values.  In  this 
respect,  the  cambered  and  symmetrical  sections  are  very  similar. 

Effect  of  Trailing-Edge  Reflex 

The  23010-1.58  airfoil  was  tested  with  its  trailing-edge  tab 
bent  upward  (reflexed)  at  an  angle  of  3  degrees.  The  and 
CM  traces  were  found  to  show  no  appreciable  differences  from 
those  for  the  airfoil  with  a  straight  trailing  edge,  aside  from 
the  small  overall  shift  in  CM-  This  agrees  with  the  steady- 
state  test  data.  It  is  therefore  concluded  that  the  favorable 
effects  noted  for  the  reflex  feature  on  cambered  blades  are  not 
associated  with  changes  in  stall  behavior. 

SUMMARY  OF  KEY  RESULTS  FROM  THE  PITCH  OSCILLATION  TESTS 

Figure  15  shows  the  highest  CN  values  attained  in  pitch  oscil¬ 
lation  for  the  Vertol  23010-1.58  airfoil.  The  Mach  0.4  and  0.6 
lines  are  significant  boundaries,  since  they  were  derived  from 
data  groups  where  changes  to  a0  did  not  significantly  increase 
the  maximum  CN .  In  the  case  of  the  Mach  0.2  line,  it  was 
apparent  that  higher  aQ's  than  those  tested  here  would  have 
given  higher  maximum  CN's  for  k's  above  0.124;  therefore,  that 
line  does  not  represent  an  actual  limit.  This  is  indicated  by 
the  change  in  slope  of  the  0.2  Mach  line  at  k  =  0.525,  where 
the  a o  =  25-degree  data  begin  to  dominate  the  aD  -  17.5-degree 
data  used  for  the  rest  of  the  M  =  0.2  line. 

The  frequency  trends  clearly  show  the  influence  of  dynamic 
effects  in  increasing  the  maximum  .  The  strong  effect  of 
Mach  number  in  suppressing  attainable  Cn  is  also  clear. 

Figure  15  shows  similar  curves  for  the  NACA  0012  (modified) 
airfoil.  The  same  frequency  trends  are  apparent,  but  at  Mach 
0 . 4  the  CN  1  imit  is  substantially  lower  than  for  the  cambered 
section . 

Damping  characteristics  of  the  Vertol  23010-1.58  airfoil  are 
presented  in  Figures  16  and  17  as  a  ratio  to  theoretical 
damping.  The  theoretical  damping  equations  are  derived  in 
Appendix  I  and  are  based  on  incompressible  nonsteady  aero¬ 
dynamic  theory.  In  Figure  16,  the  reduced  frequency  is  held 
constant  at  k  =  0.24.  Increasing  the  Mach  number  causes  the 
damping  parameter  to  cross  into  negative  values  at  lower  ac's, 
until  compressibility  alters  the  damping  characteristics  of  the 
airfoil  sufficiently  to  prevent  a  net  negative  damping  for  the 
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cycle.  At  Mach  0.6  the  lowest  value  of  negative  damping  is 
60  percent  of  the  theoretical  unstalled  value. 

In  Figure  17,  the  dimensional  frequency  is  held  constant  in¬ 
stead  of  the  reduced  frequency  as  in  Figure  16.  The  data  shown 
are  for  f  =  96  Hertz  (6/rev) ,  corresponding  to  a  typical  rotor 
blade's  torsional  natural  frequency.  This  plot  shows  that 
negative  damping  exists  only  for  a  very  narrow  range  of  Mach 
numbers  for  this  amplitude  of  oscillation;  since  the  Mach 
number  on  a  rotor  blade  varies  widely  along  the  blade's 
length,  the  contribution  to  the  total  damping  of  the  blade 
will  normally  come  from  a  limited  part  of  the  blade. 

Because  of  the  nonlinear  character  of  stall,  the  damping 
cannot  be  expected  to  be  independent  of  amplitude.  This  is 
clearly  shown  by  the  lower  part  of  Figure  18.  In  the  aQ  range 
where  instability  is  present,  larger  amplitude  motion  leads  to 
a  more  damped  condition.  This  could  lead  to  limit-cycle 
oscillations  at  the  amplitude  which  gives  zero  damping.  The 
upper  part  of  Figure  18  shows  the  influence  of  amplitude  on 
maximum  Cm .  Excursions  beyond  static  stall  are  increased  by 
oscillation  amplitude  because  of  the  higher  pitch  rate. 

The  effect  of  camber  is  illustrated  in  Figure  19,  where  the 
frequency  is  held  constant  at  a  typical  first  torsional  mode 
value.  The  beneficial  effect  of  camber  is  apparent  at  all 
Mach  numbers,  either  in  the  form  of  a  higher  aQ  before  the 
damping  becomes  negative  or  in  the  form  of  a  lower  level  of 
negative  damping. 

DRAG  IN  PITCHING  OSCILLATION 

The  problem  of  measuring  the  instantaneous  drag  of  an  airfoil 
under  unsteady-flow  conditions  is  much  more  difficult  than  that 
of  finding  the  lift.  Time-lag  and  mixing  effects  rule  out 
using  the  wake-momentum  survey  method  except  for  getting  an 
average  value.  The  use  of  a  balance  had  been  impractical  for 
Cn  and  Cm  because  of  the  complexity  of  the  dynamic  response; 
the  same  argument  would  apply  for  Cp.  It  would  be  conceivable 
to  measure  the  form  drag  from  pressure  measurements , i . e . ,  by  the 
same  technique  as  for  instantaneous  Cn •  However,  many  more 
transducers  would  be  required  because  of  the  sensitivity  of 
Cd  to  details  of  the  leading-edge  pressure.  Furthermore,  to 
obtain  Cn  and  Cc  (chord  force)  at  the  same  time,  it  would  be 
necessary  to  measure  AP's  with  respect  to  a  known  reference 
value,  instead  of  merely  taking  AP's  across  the  airfoil. 

Finally,  the  skin- f riction  component  would  not  register  in  this 
system. 

Therefore,  it  was  decided  to  measure  only  average  drag  values 
by  using  a  slowly  traversing  pitot-static  probe  to  survey  the 
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moment  defect  in  the  wake.  Drag  data  were  obtained  for  the 
Vertol  23010-1.58  airfoil  only.  All  drag  data  are  summarized 
in  Figure  20. 

The  dashed  line  shows  Cp  under  steady-a  conditions.  The  dotted 
lines  are  extrapolations  of  the  low-a  (fully  attached  flow) 
data  representing  hypothetical  Cp's  for  conditions  where  dyna¬ 
mic  effects  have  postponed  stall.  The  results  of  the  oscilla¬ 
tory  drag  measurements  show  that  the  average  Cp  goes  up  with 
frequency  between  16  Hertz  and  48  Hertz  at  both  Mach  numbers. 

Quasi-steady,  average  values  of  drag  have  been  computed  from 
the  steady  data  for  comparisons,  by  using  the  following 
equation : 


=  77  ^  CD  (ot)  da'  (8) 

where  a  is  the  instantaneous  angle  of  attack.  Two  different 
values  of  Cp  were  calculated  at  the  higher  mean  angle  of 
attack.  In  the  first  case,  the  steady-state  drag  values  were 
used  directly.  In  the  second,  the  drag  was  assumed  to  follow 
the  dotted  line  in  Figure  2C  when  da/dt  was  positive,  and  the 
steady-state  solid  line  when  da/dt  was  negative.  In  this  way 
the  effect  of  stall  delay  could  be  estimated.  At  Mach  0.4, 
it  is  seen  that  the  second  method  gives  a  more  realistic  value, 
since  the  pitching  oscillation  at  this  Mach  number  does  actually 
delay  stall  until  a  stops  increasing.  At  Mach  0.6,  the  second 
method  underpredicts  Cp.  This  is  reasonable,  since  the  simil¬ 
arity  of  the  dynamic  C^  and  Cm  traces  to  the  static  values 
implies  that  Cp  would  be  unchanged  as  well. 

Figure  20  also  shows  typical  shapes  of  total  pressure  profiles. 
It  should  be  noted  that  under  seme  conditions,  dynamic  effects 
can  actually  produce  negative  drag,  as  in  the  case  of  a  bird 
in  flight.  Such  an  effect  is  visible  in  the  center  of  the  16- 
Hertz  trace.  Further  study  of  such  phenomena  is  beyond  the 
scope  of  this  report. 

FORCED  VERTICAL  TRANSLATION 

Tests  were  run  at  frequencies  corresponding  to  16  Hertz  (1/rev) 
and  32  Hertz  (2/rev) .  Mach  numbers  from  0.2  to  0.6,  geometric 
angles  of  attack  from  0  to  20  degrees,  and  amplitudes  (Ah's) 
from  0.31  to  0.62  semichord. 

The  Cn  and  CM  traces  will  be  presented  first,  followed  by  a 
discussion  on  damping  and  the  maximum  Cjj's  attained  during  a 
cycle. 
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Effect  of  Angle  of  Attack  and  Frequency 


Figure  21  shows  the  effect  of  m0  (geometric  a,  relative  wind 
due  to  airfoil  motion  not  included)  on  and  CM  for  vertical 
translation.  Fixed  parameters  are:  Ah  =  0.30,  k  =  0.067 
(1/rev),  and  M  -  0.4.  The  data  at  ac  =  5.02  degrees  and  9.82 
degrees  show  no  sign  of  stall.  At  a  -  12.25  degrees,  and  Cm 

both  show  stall  at  the  middle  of  the  downstroke.  Stall  is 
caused  by  the  increase  in  the  effective  a  due  to  the  vector 
addition  of  plunge  rate  with  tunnel  velocity.  As  was  the  case 
for  pitching  motions,  C^  stalls  later  than  CM. 

For  plunging  motion,  the  net  work  performed  by  the  airfoil  over 
a  cycle  is  given  by 


W  =  <^CN  dh.  (9) 

This  is  proportional  to  the  area  enclosed  in  the  versus  a 
trace,  and  is  positive  (stable)  for  a  CCW  circuit. 

In  the  a 0  -  12.25-degree  case,  only  a  small  area  of  negative 
damping  appears  at  the  left  of  the  CN  trace.  At  the  next  an  , 
14.65  degrees,  the  airfoil  is  at  its  static-stal L  angle.  On 
the  upstroke,  when  dh/dt  is  positive  and  the  effective  a  is 
reduced,  the  airfoil  remains  unstalled.  At  maximum  positive  h 
and  during  the  negative  dh/dt  part  of  the  cycle,  the  airfoil  is 
fully  stalled.  This  results  in  negative  damping.  The  CM  trace 
also  shows  the  effects  of  stall.  However,  the  total  pitching 
moment  change  is  only  -0.12,  the  value  for  static  stall,  and 
no  overshoot  in  pitching  moment  is  present. 

The  traces  shown  on  Figure  22  are  at  twice  the  frequency  (2/rev) 
of  the  plots  in  Figure  21.  The  area  in  the  CN  loops  at  low 
mean  angle  of  attack  has  increased,  as  nonsteady  potential- 
flow  theory  predicts.  At  a0  =  12.46  degrees,  and  Cm  both 
show  stall  effects  beginning  at  h  =  0  on  the  downstroke.  Nev¬ 
ertheless,  negative  damping  does  not  result  because  the  air¬ 
foil  remains  stalled  until  the  middle  of  the  upstroke.  The 
stall  at  a „  =  14.88  degrees  occurs  at  the  en i  of  the  first 
quarter  of  the  cycle,  but  the  slow  variatior  of  with  time 
gives  attached  flow  only  during  the  last  ha  i  of  the  upstroke. 
There  is  a  substantial  reduction  in  damping  compared  to  the 
low-a  case,  but  it  remains  positive.  The  large  Cm's  shown 
here  may  be  significant  for  rotor  dynamics,  through  effects 
on  torsional  oscillation. 

Effect  of  Mach  Number 


Some  typical  Cm  and  Cm  traces  at  f  =  33  Hertz  (2/rev)  are  shown 
in  Figures  23  and  24  at  Mach  number  0.2  and  0.6,  respectively. 
There  is  no  evidence  of  stall  up  to  the  highest  a0  tested  at 
Mach  0.2.  A  systematic  change  is  evident  in  the  Cm  traces 
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for  M  =  0.2.  These  traces  lie  on  lines  sloping  downward  to 
the  right  and  enclose  little  or  no  area.  This  is  the  effect 
of  the  noncirculatory  loading  caused  by  plunge  acceleration. 

The  slope  of  the  line  agrees  very  well  with  the  potential 
theory  value  (dC  /dh  =  0.048  at  this  frequency).  The  tilt  was 
hardly  noticeable  at  M  =  0.4  because  it  is  proportional  to  the 
square  of  the  reduced  frequency.  At  M  =  0.4,  dCM/dh  would 
have  been  only  -0.012. 

At  M  =  0.6,  there  is  a  substantial  amount  of  load  rearrange¬ 
ment,  as  evidenced  by  the  moment  hysteresis  loops.  At  this 
Mach  number,  the  effect  of  increasing  a  is  to  flatten  the 
potential-flow  ellipse  on  the  CN  versus  h  plot.  This  severely 
limits  the  variation  of  lift  with  h.  No  figure-eight  features 
of  appreciable  magnitude  appear  in  either  the  CN  or  the  Cm 
traces;  this  can  be  ascribed  to  the  absence  of  sharp  breaks  in 
the  static  data,  as  in  the  case  of  pitching  motion. 


Effect  of  Amplitude 


Figure  25  shows  the  effect  of  a  50-percent  increase  in  trans¬ 
lation  amplitude  on  the  C  and  C  versus  a  traces.  From  com¬ 
parison  of  these  loops  with  the  third  set  of  data  in  Figure 
21,  it  is  clear  that  the  general  characteristics  of  the  data 
are  preserved  when  the  amplitude  is  increased. 


Effect  of  Airfoil  Symmetry 


Figure  26  shows  Cn  and  Cm  traces  for  the  NACA  0012  (modified) 
airfoil  in  plunging  motion  at  32  Hertz  (h  =  0.24),  Mach  0.2, 
and  Ah  =  0.31.  At  aQ  =  10  degrees  (well  below  stall),  the 
traces  are  practically  identical  to  those  for  the  23010-1.58 
airfoil.  At  12.5  and  15  degrees,  however,  very  substantial 
stall  effects  appear;  the  C  loop  becomes  somewhat  irregular, 
its  center  sagging  to  a  level  of  1.0  or  less,  and  a  large  drop 
in  C  appears  at  the  end  of  the  downstroke.  This  is  in  strong 
contrast  to  the  behavior  of  the  cambered  airfoil  under  the 
same  conditions.  Data  at  M  =  0.4  are  not  shown,  since  var¬ 
iance  from  the  behavior  of  the  23010-1.58  airfoil  was  less 
striking,  although  stall  effects  were  more  prevalent.  At 
Mach  0.6,  no  noticeable  differences  were  found  between  the 
airfoils . 
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SUMMARY  OF  KEY  RESULTS  OF  VERTICAL  TRANSLATION  TESTING 


Maximum  Normal  Force  Data 


Figure  27  shows  the  maximum  Cjg '  s  obtained  during  tests  in  the 
plunging  motion  mode.  The  values  reached  were  very  insensi¬ 
tive  to  the  frequency  and  amplitude  of  motion,  though  the 
maxima  were  generally  reached  at  conditions  corresponsing  to 
the  largest  rate  of  change  of  instantaneous  i  (with  respect  to 
the  relative  wind).  At  Mach  0.2,  no  maximum  was  reached  for 
the  23010-1.58  section  for  the  i0's  tested.  The  value  shown 
is  extrapolated. 

Both  airfoils  show  a  increase  over  the  steady-flow  value, 
ana  the  effect  of  larger  Mach  number  is,  as  usual,  to  suppress 
the  attainable.  The  higher  stall  angle  of  the  23010-1.58 
airfoil  is  apparent  in  its  consistent  increment  of  0.2  in 
over  the  symmetrical  section. 

Damping  in  Plunge 

The  tests  reported  here  are  the  first  known  to  the  authors  in 
which  negative  average  damping  was  observed  in  translation 
motion.  These  tests  were  conducted  at  a  much  larger  Ah  and 
at  higher  Mach  numbers  than  any  previous  tests.  This  is  pro¬ 
bably  why  neither  Halfman  *  nor  Rainey?  measured  negative 
damping . 

Figure  28  shows  damping  in  plunge  for  the  Vertol  23010-1.58 
airfoil  section  as  a  function  of  geometric  a Q  for  M  =  0.4  and 
M  =  '0 .  G .  For  angles  well  below  the  static-stall  level,  damping 
at  M  =  0.4  remains  near  the  theoretical  value.  This  is  con¬ 
sistent  with  lift  variation  at  the  rate  corresponding  to  a 
stall-free,  lift-curve  slope.  The  Loginning  of  the  loss  of 
damping  is  apparent  at  12.3  degrees;  above  15  degrees,  zero  or 
negative  damping  prevails.  The  largest  negative  value  occurs 
when  i0  is  just  at  the  stall  point  and  the  plunging  motion 
varies  the  aerodynamic  a  just  enough  to  cross  the  lift-curve 
discontinuity.  (The  CN  and  CM  traces  for  the  cycle  are  shown 
in  Figure  21 .  ) 

At  Mach  0.6,  the  loss  of  damping  occurs  at  a  lower  uD,  but  no 
large  negative  values  were  observed.  This  is  consistent  with 
a  CN  versus  a  relationship  in  which  lift  has  practically 
ceased  to  vary.  This  agrees  with  the  steady-flow  Cn  data  and 
the  Cn  and  Cm  versus  a  curves  of  Figure  52.  No  stall  was 
observed  for  this  profile  at  Mach  0.2.  Therefore,  the  damping 
was  always  observed  to  be  near  the  theoretical  level  and  is 
not  shown. 


Figure  29  shows  similar  data  for  the  symmetrical  airfoil.  The 
damping  drops  off  at  a  lower  a0  than  for  the  cambered  profile 
at  Mach  0.4,  but  near  the  same  poin..  at  Mach  0.6.  The  maximum 
negative  value  reached  was  only  half  that  of  the  NACA  23010-1.58 
airfoil.  Presumably,  this  is  a  consequence  of  the  smaller 
discontinuity  in  the  C„,  versus  a  curve. 

Both  the  symmetrical  and  the  cambered  airfoils  show  progres¬ 
sively  less  severe  negative  damping  for  larger  h  excursions. 

This  implies  the  possibility  of  limit-cycle  bending  oscilla¬ 
tions  of  rotor  blades. 

DYNAMIC  EFFECTS  ON  CHORDWISE  LOADING 

Considerable  insight  into  the  processes  governing  stall,  C 
and  CM  breaks,  and  negative  damping  can  be  gained  by  examining 
ACp  time  histories  and  chordwise  load  distributions.  To  pro¬ 
vide  a  background  for  comparison  with  dynamic  results,  and  to 
explain  the  existence  of  separate  branches  on  the  Cjj  and  Cm 
versus  a  curves,  chordwise  loadings  for  steady-flow  conditions 
will  be  examined  first. 

Chordwise  Loadings  Under  Static  Conditions 

Figure  30  shows  how  Mach  number  determines  whether  c  iscontinu- 
ities  will  be  present  on  the  moment  and  normal  force  curves. 
Pressure  distributions  are  shown  for  angles  of  attack  just 
below  and  just  above  the  Cm  or  Cm  break  at  Mach  0.4  and  at 
Mach  0.6.  At  Mach  0.4,  the  loading  is  seen  to  be  of  the 
classical  potential-flow  type  before  stall.  The  large  negative 
d(ACp)/dx  implies  that  the  upper  surface  pressure  must  rise 
rapidly  just  aft  of  the  leading  edge.  The  result  of  a  slight 
increase  in  a  is  to  aggravate  this  adverse  pressure  gradient, 
precipitating  stall.  The  redistribution  of  load  is  drastic: 

Cn  drops  by  20  percent  and  CM  falls  to  about  0.1.  The  flat¬ 
ness  of  the  ACp  curve  is  symptomatic  of  a  large  region  of 
separated  flow,  filled  with  turbulent  eddies. 

Compare  those  results  with  the  Mach  0.6  loadings:  a  completely 
different  phenomenon  operates  in  this  case.  The  pressure  dis¬ 
tributions  for  both  angles  of  attack  imply  a  favorable  pres¬ 
sure  gradient  forward  of  10-percent  chord.  This  is  due  to 
transonic  effects.  Supersonic  expansion  around  the  leading 
edge  causes  increasing  local  air  velocity  and  falling  pres¬ 
sure  until  a  shock  is  reached  near  15-percent  chord. 

When  a  increases,  the  flow  is  able  to  adjust  to  the  change  by 
a  small  readjustment  of  the  expansion/shock  system.  Slightly 
weakened  pressure  recovery  results,  leading  to  small  increases 
in  negative  pitch  moment  and  in  C^.  The  Mach  0.6  flow  is  thus 
free  of  discontinuities,  though  changes  in  slope  in  dCfj/da 
and  dC,./da  do  occur. 
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Loadings  During  Pitching  Oscillations 


Effect  of  Frequency 

Figure  31  shows  pressure,  normal  force,  and  pitching  moment 
variations  with  time  for  the  Vertol  23010-1.58  airfoil 
oscillating  at  16  Hertz.  Even  at  this  comparatively  low 
frequency,  substantial  postponement  of  stall  is  evident 
on  the  C..  and  C,.  curves.  Note  also  that  the  AC*,  nearest 
the  leading  edge  (0.01c)  remains  substantially  greater 
than  the  AC  further  down  the  chord  throughout  the  period 
of  stall. 

Figure  32  shows  a  sequence  of  chordwise  loadings  at  20- 
degree  intervals  in  the  cycle.  Theoretical  steady-flow 
loading  profiles  are  shown  for  comparison.  Note  that 
the  measured  loadings  match  the  theoretical  shapes  with 
great  precision  when  stall  is  not  present. 

Figures  33  and  34  show  data  similar  to  those  of  Figures  31 
and  32,  except  that  the  frequency  here  is  96  Hertz.  The 
dynamic  influence  on  Cn  and  Cm  is  dramatic.  The  Cm  trace 
now  forms  a  passable  sinusoid,  not  very  different  from 
the  quasi-steady,  potential-flow  line.  However,  the 
presence  of  stall  is  evident  on  the  CM  curves.  The  chord- 
wise  loadings  show  a  slight  but  noticeable  difference  from 
the  steady-flow  theoretical  shapes  over  the  increasing-u 
range.  This  is  caused  by  effective  additional  camber  due 
to  pitch  rate. 

Effect  of  Mach  Number 

It  has  been  noted  that  particularly  violent  fluctuations 
of  force  and  moment  occurred  when  dynamic  effects  were 
strong  enough,  because  of  low  Mach  number,  to  delay  stall. 
An  exceptionally  large  CN  was  reached  in  this  case. 

Figure  35  shows  the  loading  variation  on  the  Vertol  23010- 
1.58  airfoil  as  the  flow  breaks  down  at  a  =  22.5  degrees 
(16  Hertz,  Mach  0.2).  The  differential  pressure  across 
the  aft  portion  of  the  blade  is  three  times  larger  than 
for  steady-flow  stall. 

The  effect  of  Mach  number  on  high-frequency  oscillatory 
loadings  is  shown  by  Figures  36  through  39,  which  form 
a  complete  sequence  when  considered  with  Figures  31  and 
32.  At  Mach  0.2,  the  pressure  profiles  shown  in 
Figure  37  are  completely  free  of  visible  signs  of  stall. 

At  every  a,  ACp  falls  steadily  toward  the  aft  portion  of 
the  airfoil. 

Oscillations  at  96  Hertz  from  Mach  0.4  have  already  been 
discussed.  In  that  case,  dynamic  effects  were  not  able 
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to  eliminate  substantial  stalling.  At  Mach  0.6,  the 
situation  is  again  different.  Good  pressure  recovery 
appears  to  prevail  under  all  conditions.  The  influence 
of  nose-dvwn  pitch  rate  as  effective  inverse  camber  is 
discernible  in  the  loading  profiles  corresponding  to 
nose-down  motion.  Then,  the  expansion/shock  effect  is 
exaggerated,  and  almost  complete  unloading  of  the  aft 
part  of  the  section  results  from  this  high  Mach  number. 
Steady-flow  tests  often  show  such  loadings.  (See 
Reference  8  for  examples  of  this.) 


Effect  of  Airfoil  Symmetry 


Figures  40  and  41  show  time  histories  and  chordwise 
loadings  for  the  NACA  0012  (modified)  airfoil  oscil¬ 
lating  at  80  Hertz  in  a  Mach  0.4  flow.  The  most  striking 
difference  in  the  results  shown  for  this  cambered  sec¬ 
tion  is  the  complete  loss  of  the  leading-edge  pressure 
peak  when  stall  occurs.  This  is  characteristic  of  the 
symmetrical  airfoil  under  all  conditions  tested.  It  is 
partly  responsible  for  the  consistently  larger  negative 

stalled  C, . 1  s . 

M 

Note  also  how  a  ACp  wave  appears  to  travel  down  the  air¬ 
foil  in  the  period  Er^m  80  to  200  degrees  of  the  cycle. 
This  wave  could  be  the  result  of  the  shedding  of  a  strong 
vortex  from  the  leading  edge,  which,  in  turn,  is  a  con¬ 
sequence  of  the  loss  of  pressure  peak  already  discussed. 
This  is  the  other  contributor  to  the  increased  negative 
CM  values  observed  for  the  symmetrical  section. 


Loadings  in  Translation  Motion 


Figures  42  and  43  show  typical  effects  of  plunging  oscillation 
on  the  cambered  airfoil  loadings,  while  Figures  44  and  45  give 
corresponding  data  for  the  symmetrical  airfoil.  The  angle  of 
attack  selected  for  the  NACA  0012  plot  is  deliberately  lower, 
since  approximately  equal  margins  from  stall  were  desired. 


The  results  in  both  cases  are  consistent  with  those  observed 
for  pitching  oscilla cions .  Note  especially  the  wider  Cm 
excursion  of  the  0012,  and  also  its  more  serious  loss  of 
leading-edge  suction. 


COMPARISON  WITH  THEORETICAL  RESULTS 


Experimental  force  and  moment  derivatives  for  both  the  symme¬ 
trical  and  the  cambered  airfoils  are  compared  with  theore¬ 
tically  predicted  values  in  Figure  46.  The  theoretical  data 
and  the  wind  tunnel  wall  corrections  were  calculated  with  a 
computer  program  which  was  supplied  by  NASA^.  The  theory  was 
based  on  thin  airfoil  representation  with  the  effects  of 
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compressibility  and  the  tunnel  walls  (method  of  images)  in¬ 
cluded.  Test  data  are  presented  only  for  small  a0's,  since 
the  viscous  and  transonic  flow  effects  which  appear  at  high 
values  of  aQ  cannot  be  predicted  analytically. 

Pitching  Motion 

The  experimental  values  of  the  dynamic-lift  curve  slope  are 
lower  than  the  theoretical  ones  in  all  cases.  This  agrees  with 
the  universally  noted  reduction  of  the  static-lift  curve  slope 
caused  by  viscosity.  The  trends  to  higher  values  of  | Cm ( / | a | 
with  increasing  Mach  number  agree  well  with  the  predictions. 

The  Cm  phase  agrees  well  at  Mach  0.2,  but  at  0.4  and  G.6  the 
experimental  results  lag  by  as  much  as  18  degrees.  No  expla¬ 
nation  for  this  discrepancy  is  now  available.  More  of  such 
lag  is  apparent  in  the  case  of  the  Vertol  23010-1.58  than  for 
the  NACA  0012  (modified) .  This  is  caused  by  separation  effects 
on  the  lower  surface  near  the  leading  edge. 

Pitching  moment  derivatives  also  show  generally  good  agreement, 
but  again  they  show  some  reduction  in  magnitude.  The  same 
thing  was  noted  by  Halfman^  and  is  attributed  to  viscosity. 

The  trend  with  Mach  number  is  identical  for  both  experimental 
and  theoretical  data.  The  Cm  phase  of  the  symmetrical  airfoil 
moment  agrees  with  theory  within  7  degrees.  For  the  cambered 
airfoil,  the  difference  between  theory  and  test  is  20  degrees, 
at  low  reduced  frequency.  This  is  caused  by  a  slight  rearward 
shift  of  the  aerodynamic  center  from  the  quarter-chord  point 
f^r  the  additional  loading  due  to  oscillations.  The  20-degree 
difference  in  phase  represents  a  shift  of  only  0.8-percent 
chord,  owing  to  the  extreme  sensitivity  of  the  pitching  moment 
phase  to  the  location  of  the  pitch  reference  center. 

Vertical  Translation 


The  test  data  for  Cjg  and  CM  in  translation  are  compared  with 
theory  in  Figure  47.  The  agreement  of  the  Cn  and  Cm  slopes 
with  oscillatory  a  is  very  good  except  at  very  low  values  of  k 
for  the  cambered  airfoil,  where  the  Cm  magnitude  response  is 
aga_n  much  larger  than  theory  predicts.  This  can  be  expected, 
since  the  behavior  of  the  cambered  airfoil  at  k's  of  0.03  will 
be  essentially  quasi-steady,  and  the  aerodynamic  center  of  the 
cambered  airfoil  is  not  on  the  quarter-chord. 

There  is  a  20-degree  difference  in  the  oscillating  CN  phase 
between  theory  and  test  for  the  symmetrical  airfoil.  There  is 
also  a  15-degree  difference  between  the  C^j  phases  of  the 
symmetrical  and  cambered  airfoils.  No  such  difference  was 
noted  for  Cm;  no  explanation  for  this  is  known  at  the  present 
time.  The  difference  between  theory  and  test  for  the  moment 
phase  can  be  explained  by  an  offset  of  the  center  of  pressure 
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from  the  quarter -chord.  Again,  the  Cm  phase  is  very  sensitive, 
if  moments  are  taken  about  an  axis  near  tho  aerodynamic  center. 

COMPARISON  WITH  PREVIOUS  RESULTS 

United  Aircraft  Corporation  Oscillation  Data 

Average  pitch  damping  per  cycle  is  compared  with  Carta's 
results  in  Figure  48.  The  range  of  a0  for  negative  damping 
agrees  well,  but  the  tests  reported  here  indicate  a  lower 
minimum  value  than  do  the  UAC  experiments. 

The  test  program  reported  here  had  the  benefit  of  automated 
data  reduction  methods.  Therefore,  all  of  the  results  shown 
are  averaged  from  at  least  five  consecutive  cycles.  In  the 
case  of  the  symmetrical  airfoil  data,  ten  cycles  were  used. 

The  UAC  results  were  manually  processed,  so  only  one  cycle  per 
test  point  could  be  covered.  Therefore,  differences  are  to 
be  expected. 

CH-47A  Dynamic  Airloads  Data 

Cm  data  from  the  dynamic  airloads  tests10  are  compared  with 
symmetrical  airfoil  data  from  this  report  in  Figure  49.  The 
tunnel  test  point  was  chosen  to  represent  the  Mach  and  a 
conditions  on  the  retreating  side  of  the  rotor  disc  (M  =  0.4, 
ae  =  10  degrees  +  5  sin  0)  and  the  pitch  reference  position 
has  been  shifted  to  agree  with  the  stall  point  of  the  flight 
test  data.  The  agreement  between  test  and  theory  is  remark¬ 
able,  both  in  the  magnitude  of  the  Cm  change  with  stall  and  in 
the  duration  of  stall.  Both  the  wind  tunnel  and  the  flight 
test  data  show  a  partial  reattachment  at  the  270-degree  azi¬ 
muth  position  followed  by  a  further  stall.  The  flight  test  data 
show  a  6  to  7/rev  oscillation  on  the  advancing  side,  which  is 
probably  elastic  response  in  the  first  torsion  mode.  The 
nose-down  pitching  moment  of  the  flight  test  data  in  the  60- 
to  180-degree  azimuth  range  may  be  due  to  the  higher  Mach 
number  for  those  data. 

TIME  HISTORIES  OF  STARTS  AND  STOPS  FOR  PITCHING  OSCILLATION 

A  complete  time  history  of  forced  pitch  oscillation  from 
beginning  to  end  was  recorded  for  ten  data  points.  Since  the 
data  from  these  tests  span  many  records  and  many  cycles  of 
oscillations,  they  are  not  reducible  by  the  computer  programs 
developed  for  the  airfoil  in  steady  oscillation;  therefore, 
they  are  not  included  in  the  data  listings. 

Time  history  plots  from  the  analog  traces  have  been  scanned  for 
unusual  behavior.  Nothing  worthy  of  note  has  been  found. 
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PITCH  OSCILLATION  OF  AN  AIRFOIL  SYSTEM  WITH  A  TORSIONAL  DEGREE 
OF  FREEDOM 


Because  of  the  widely  varying  aerodynami c  parameters,  such  as 
the  reduced  frequency  and  Mach  number,  along  a  rotor  blade  in 
forward  flight,  the  stall-induced  torsional  oscillations  of  a 
whole  blade  cannot  be  simulated  in  a  wind  tunnel.  Another  com¬ 
plication  is  the  fact  that  the  blale  is  not  executing  simple 
sinusoidal  pitching  motion  at  the  conditions  of  greatest  interest. 
Rather,  it  is  producing  a  combination  of  the  l/'rev  basic  angle- 
of-attack  change  (discussed  in  the  blade  motion  studies  section) , 
with  varying  amounts  of  higher  frequency  elastic  notions  super¬ 
imposed.  These  elastic  pitching  and  translation  deflections 
modify  the  angle  of  attack.  Of  these  elastic  motions  the  res¬ 
ponse  of  the  first  torsional  mode  is  of  most  interest,  since  it 
can  be  forced  tc  large  deflections  by  stall  effects. 


Therefore,  the  approach  taken  in  this  series  of  tests  was  to 
simulate  a  representative  blade  section  by  spring-mounting  the 
airfoil,  which  was  tuned  to  4  to  6/rev  (64  to  96-Hertz  model 
scale) ,  and  by  driving  the  angle  of  attack  at  1/rev  (16-Hertz 
model  scale)  with  the  system  used  for  forced  pitching  tests. 


Oscillations  in  the  first  torsional  mode  for  a  rotor  blade  are 
governed  by  a  second-order  differential  equation  of  the  form 


a"  +  a 


l/2p ( 2b)  (area) 


CM  * a  '  a 


.)  . 


(10) 


The  conditions  required  for  aerodynamic  similarity  are: 


1.  A  correct  Mach  number  to  simulate  the  aerodynamic 
environment  of  the  retreating  blade. 

2.  A  Reynolds  number  in  the  two-  to  six-million  range 

3.  A  reduced  frequency  corresponding  to  rotor  1/rev  for 
the  background  pitch  motion  contribution  to  a,  a',  etc. 

4.  The  elastic  torsional  frequency  must  be  four  to  six 
times  the  background  frequency. 


The  condition  for  dynamic  similarity  requires  that  the  coeffi¬ 
cient  of  CM  on  the  right-hand  side  of  the  equation  for  the 
rotor  blade  motion  match  the  corresponding  coefficient  in  the 
equation  for  airfoil  system  motion.  That  coefficient  is  called 
the  torsional  Locke  number  and  is  denoted  here  by 

yt  can  also  be  written  as 


2ob4 

f  v  1 

2  =  20b11 

1 

T 

IZ F 

I  _ 

~TT 

T 

l  T 

T 

T 

9 


35 


The  factor  kT,  the  reduced  torsional  natural  frequency,  is 
dictated  by  aerodynamic  similarity  and  can  be  simulated  easily 
in  a  tunnel.  It  was  determined  by  model  strength  requirements, 
and  the  resulting  value  of  2pb^  (span) /It  was  low  by  a  factor 
of  10.  The  result  of  this  difference  is  a  reduction  to  the 
model's  displacement  to  10  percent  of  the  full-scale  level  due 
to  a  given  aerodynamic  force. 


While  the  model  tests  will  not  provide  full-scale  deflections, 
they  can  provide  the  following  important  results: 


1 .  The  tests  can  predict  the  onset  of  stall-induced 
oscillations,  because  these  depend  on  the  1/rev  a 
variation . 


2.  The  tests  can  provide  data  with  which  analytical 

formulations  of  the  effects  of  stall  on  rotor  blades 
can  be  checked  in  a  dynamic  system  described  by  the 
torsional  equation. 


The  pitch  response  data  from  the  tuned  system  tests  are  pre¬ 
sented  in  Figures  50  and  51.  Only  the  Vertol  23010-1.58 
airfoil  was  tested  in  this  manner  for  the  test  conditions 
listed  in  Tables  IX  and  X  in  Volume  II. 


Effect  of  Mean  Angle  of  Attack  and  Mach  Number 


The  amplitude  of  the  maximum  elastic  pitch  response  is  directly 
related  to  the  angle  for  static  stall  as  shown  in  Figure  50  for 
Mach  numbers  0.4  and  0.6.  The  maximum  higher  harmonic  6a, 
rather  than  the  fourth,  fifth,  or  sixth  harmonic  of  the  tuned 
system,  has  been  plotted  because  the  peak  response  frequency  of 
the  system  is  affected  by  damping  and  by  the  slight  variations 
in  the  1/rev  driving  frequency.  All  pitch  response  data  have 
been  multiplied  by  kT  ,  which  is  the  dominant  term  in  Locke 
number  (b^/lT  is  constant  for  the  model) .  The  elastic  response 
is  decreased  by  changing  the  Mach  number  from  0.4  to  0.6.  This 
agrees  with  the  conclusions  drawn  from  forced  pitch  oscillation 
tests . 


The  higher  harmonic  response  for  and  Cm  is  in  general  agree¬ 
ment  with  the  aQ  response.  At  mean  angles  of  attack  where  the 
motion  does  not  experience  stall  effects,  the  higher  harmonic 
pitch  response  rapidly  approaches  the  background  levels  measured 
in  the  forced  pitch  oscillation  testing. 
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Effect  of  System  Resonant  Frequency 

Figure  50  shows  the  effect  of  changing  the  system  resonant 
frequency  from  64  to  99  Hertz  (4  to  6/rev) .  The  increase  in 
response  is  spectacular  for  the  4/rev  system.  This  shows  that 
the  response  to  stall  is  highly  dependent  on  the  stiffness  of 
the  system. 

Effect  of  Driving  Amplitude 

Two  pitch-driving  amplitudes  were  tested  (a  =  +5  degrees  and 
i  =  +  7.5  degrees).  The  effect  on  pitch  response  is  shown  in 
Figure  51  for  M  =  0.4  at  a  tuned  frequency  of  f^  =  80  Hertz. 
The  2.5-degree  a Q  shift  between  the  two  curves  indicates  a 
primary  dependence  of  the  pitch  response  on  the  stall  behavior 
of  the  motion  and  not  on  the  driving  amplitude. 

Effect  cf  Locke  Number 


A  few  test  points  were  obtained  at  50-percent  total  pressure  in 
the  test  section.  This  test  has  the  combined  effects  of 
reducing  the  Reynolds  number  and  the  Locke  number,  y^.  Since 
the  forced  pitch  oscillation  data  did  not  show  differences  with 
Reynolds  number,  the  primary  effect  of  the  lower  pressure  was 
to  change  yx-  In  all  cases,  the  elastic  a  response  was 
decreased  with  a  decrease  in  y  . 

Effect  of  External  Damping 

Two  levels  of  external  damping  were  used  during  the  tests.  The 
effect  of  external  damping  is  shown  in  Figure  51  for  the  aQ  for 
maximum  response.  The  maximum  elastic  response  was  decreased 
by  a  factor  of  three  with  the  application  of  full  damping  (c  = 
0.10)  . 
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CONCLUSIONS 


1.  Pitching  motions  at  the  amplitudes  and  frequencies  topical 
of  blade  elements  of  helicopter  rotors  strongly  affect  the 
limits  of  aerodynamic  normal  force  due  to  stall.  The  maxi¬ 
mum  normal  force  coefficient  in  all  cases  shows  an 
increase  which  grows  with  frequency  but  which  declines  as 
the  Mach  number  approaches  0.6. 

2.  The  margin  of  maximum  normal  force  for  cambered  over 
symmetrical  airfoils  which  has  been  observed  in  the  past 
under  steady-flow  conditions  prevails  under  dynamic 
pitching  conditions  as  well. 

3.  Dynamic  effects  postpone  symptoms  of  stall  in  airfoil  pitch 
ing  moment  also.  However,  the  increase  in  angle  of  attack 
for  development  of  substantial  nose-down  pitching  moment 

is  generally  smaller  than  the  increase  in  a  for  normal 
force  break. 

4.  Pitching  oscillations  at  angles  of  attack  where  flow  separa 
tion  occurs  can  result  in  a  net  input  of  energy  from  the 
passing  airstream  to  the  airfoil  dynamic  system.  (In  the 
case  of  a  rotor  blade,  this  would  imply  increasing  ampli¬ 
tude  of  oscillation.)  This  negative  damping  diminishes 
for  larger  amplitudes  of  oscillation.  It  is  sensitive  to 
Mach  number,  but  it  disappears  for  all  practical  purposes 
on  increasing  Mach  number  to  Mach  0.6. 

5.  The  cambered  airfoil  shows  an  advantage  of  3  to  6  degrees 
in  angle  of  attack  over  the  symmetrical  airfoil  at  Mach 
numbers  where  negative  damping  can  occur.  Camber  also 
reduces  the  magnitude  of  the  pitching  moments  developed. 

6.  Oscillatory  vertical  translation  (plunge)  can  result  in 
increased  maximum  normal  force  over  steady-flow  conditions. 

7.  The  steady-flow  advantage  of  camber  in  maximum  normal  force 
coefficient  carries  over  into  dynamic  stall  in  plunge. 

8.  Plunging  motion  at  the  amplitudes  and  frequencies 
pertinent  to  helicopter  rotor  blades  produces  negative 
damping  over  narrow  ranges  of  angle  of  attack  at  Mach 
0 . 2  and  Mach  0.4. 

9.  The  effects  of  negative  damping  in  pitch  and  plunge  kre 
strongly  related  to  discontinuities  in  the  static  normal 
force  and  pitching  moment  curves.  These  discontinuities 
are  caused  by  leading-edge  stall.  At  Mach  0.6,  where 
transonic  effects  eliminate  leading-edge  stall,  damping  is 
reduced  below  the  potential  flow  level  but  remains  positive 
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RECOMMENDATIONS 


1.  The  stall  behavior  data  presented  in  this  report  should  be 
analyzed  with  a  view  to  mathematical  description.  Such  a 
description  is  a  necessary  first  step  in  the  practical 
application  of  these  data  to  rotor  dynamic  and  aerodynamic 
analyses. 

2.  Thin  airfoils  as  proposed  for  advanced-technology  helicopter 
rotors  should  be  tested  under  dynamic  conditions  to  estab¬ 
lish  their  damping  and  normal  force  characteristics. 

3.  A  program  of  analysis  should  be  undertaken  to  develop  a 
detailed  understanding  of  the  flow  processes  affecting 
dynamic  stall,  such  as  local  supersonic  flow  and  vortex 
shedding.  The  program  would  require  supplemental  experi¬ 
mental  work  in  flow  visualization  and  the  measurement  of 
absolute  pressures  to  define  local  flow  Mach  number  on  the 
airfoil  surfaces. 
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gure  9.  Effect  of  Frequency  on  Cfj  and  Tj.rae  History 
at  Stall  for  Vertol  23010-1.58  Airfoil  at 
M  =  0.4,  a 0  =  12.3°,  and  Aa  =  5°. 
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PITCHING  MOMENT  COEFFICIENT  (Cm)  NORMAL  FORCE  COEFFICIENT  <CK) 


ANGLE  OF  ATTACK  (a]  -  DEGREES 

Figure  10.  Amplitude  of  Oscillation  Effect  on  Vertol 

23010-1.58  Airfoil  at  M  =  0.4  and  f  =  32  Hertz. 
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PITCHING  MOMENT  COEFFICIENT  (CM)  NORMAL  FORCE  COEFFICIENT  <CN) 


PITCHING  MOMENT  COEFFICIENT  (CM)  NORMAL  FORCE  COEFFICIENT  <CN) 


t 


ANGLE  OF  ATTACK  (a)  -  DEGREES 


ANGLE  OF  ATTACK  (a)  -  DEGREES 


ANGLE  OF  ATTACK  ( a )  -  DEGREES 


Figure  12.  Dynamic  Cjj  and  CM  Versus  a  for  Vertol  23010-1.58 
Airfoil  at  DJ  =  0.6,  f  =  16  Hertz,  and  Aa  =  5°. 
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NOTES:  VERTOL  23010-1.58 


1.  <*0  =  17.5°  and  25°  (M  =  0.2) 

2.  a0  *  17.5°  (M  =  0.4) 

3.  a0  =  12.5°  (M  =  0.6) 

NACA  0012  (Modified) 

4.  a0  -  17.5°  and  25°  (M  =  0.2) 

5.  7.5  <  a0  <  17.5°  (M  -  0.4) 

6.  a0  =  10°  (M  =  0.6) 


REDUCED  FREQUENCY  -  k 


Figure  15.,  Maximum  Normal  Force  Attained  During  These 
Tests  for  Pitch  Oscillation  at  Aa  =  5°. 
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DAMPING/THEORETICAL  DAMPING 


NOTES :  1 . 
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a  0  -  DEGREES 


Figure  16.  Effect  of  Mach  Number  on  the  Aerodynamic  Pitch 
Damping  of  Vertol  23010-1.58  Airfoil. 
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Figure  17.  Effect  of  Mach  Number  on  the  Aerodynamic  Pitch 

Damping  of  Vertol  23010-1.58  Airfoil  at  Constant 
Frequency. 


Figure  18.  Amplitude  of  Oscillation  E 
Damping  for  Vertol  23010-1 
at  M  =  0.4  and  f  =  16  Hert 


a  =  5°  for  Vertol  23010-1.58  and 
0012  (Modified)  Airfoils. 
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TRANSLATION  POSITION  (h) 


TRANSLATION  POSITION  (h) 


Figure  23.  Dynamic  CN  and  CM  Versos  h  for  Vertol  23010-1.58 
Airfoil  at  M  =  0.2,  f  =  33  Hertz,  and  Ah  =  0.306. 
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Figure  26.  Dynamic  and  CM  Versus  h  for  NACA  0012  (Modified) 

Airfoil  at  M  =  0.2,  f  =  32  Hertz,  and  Ah  =  0.306. 
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Figure  28.  Summary  of  Damping  Data  in  Vertical  Translation 
for  Vertol  23010-1.58  Airfoil. 
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Summary  of  Damping  Data  in  Vertical  Translation 
for  NACA  0012  (Modified)  Airfoil. 
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NOTES:  1.  TP  ■  1102.4 


2  .  M  -0.4 

3.  f  -  94.34  Hz  (k  -  0. 365) 


AIRFOIL  PITCH  REFERENCE!  )  -  DEGREES 


Figure  33.  Stall  Time  History  of  CN,  CM,  and  ACp  for 
Vertol  23010-1.58  Airfoil  in  Pitch  at 
M  =  0.4,  f  =  94.3  Hertz,  and  a  =  12.5°. 
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Figure  35.  Load  Distributions  on  the  Vertol  23010-1.58 
Airfoil  During  a  High-a  Dynamic  Stall  at 
M  =  0.2,  f  =  16  Hertz,  and  Aa  =  5°. 
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NOTES 


2.  M  -  0.2 


AIRFOIL  PITCH  REFERENCE  (  )  -  DEGREES 

Figure  36.  Stall  Time  History  of  CN ,  CM,  and  aCp  for 
Vertol  23010-1.58  Airfoil  in  Pitch  at 
M  =  0.2,  f  =  96.6  Hertz,  and  ao  =  15.0°. 
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Figure  38.  Stall  Time  History  of  CN ,  CM,  and  ACp  for 
Vertol  23010-1.58  Airfoil  in  Pitch  at 
M  =  0.6,  f  =  98.0  Hertz,  and  ao  =  7.5°. 
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NOTES :  1 .  TP  -  3200 . 4 

2.  M  -  0.4 


Figure  40.  Stall  Time  History  of  CN,  Cw,  and  ACp  for 
NACA  0012  (Modified)  Airfoil  in  Pitch  at 
M  =  0.4,  f  =  89.3  Hertz,  and  aQ  =  12.2°. 
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Figure  42.  Stall  Time  History  of  ,  CM,  and  ACp  for 
Vertol  23010-1.58  Airfoil  in  Translation 


at  M  =  0.4,  f  =  30.1  Hertz,  and  a0  =  12.5°. 
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NOTES:  1.  TP  =  4041.1 

2.  M  =  0.4 

3.  f  =  30.12  Hz 

(k  =  0.116) 


a  Q  =  12.45 


4. 

5.  Afi.  =  0.47 
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Figure  43.  Stall  Time  History  of  Chordwise  Pressures  for 
Vertol  23010-1.58  Airfoil  in  Translation  at 
M  =  0.4,  f  =  30.1  Hertz,  and  aQ  =  12.5°. 
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NOTES : 


1 .  TP  *  3107.1 

2 .  M  “  0.4 


AIRFOIL  PITCH  REFERENCE (  )  -  DEGREES 


Figure  44.  Stall  Time  History  cf  ,  CM ,  and  ACD  for 

NACA  0012  (Modified)  Airfoil  in  Translation 
at  M  =  0.4,  f  =  30.3  Hertz,  and  an  =  12.2°. 
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Figure  45-  Stall  Time  History  of  Chordwise  Pressures  for 
NACA  0012  (Modified)  Airfoil  in  Translation  at 
M  =  0.4,  f  =  30.3  Hertz,  and  a0  =  12.2°. 
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Figure  48. 
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Comparison  of  Damping  with  Data  of  Reference  4 
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Fiqure  50.  Effect  of  Mach  Number  and  System  Resonant 
Frequency  on  Elastic  Pitch  Response. 
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APPENDIX  I 

DERIVATION  OF  THEORETICAL  DAMPING 


In  order  to  nondimensionalize  the  cycle  damping  parameter,  the 
corresponding  value  for  the  incompressible  case  will  be  derived 
using  the  Theodorsen1  formulation  for  normal  force  and  moment: 


L  =  +  L2  + 


(12) 


where 


is  circulatory  origin,  acting  at  quarter-chord 

L2  is  noncirculatory  origin,  acting  at  mid-chord 

L3  is  noncirculatory  origin,  acting  at  three-quarter- 
chord 
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The  equations  can  be  reduced  to  coefficient  form  by  dividing 
by  l/2pV2(2b)  (1)  with  V  assumed  to  be  a  constant. 
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Similarly,  the  pitching  moment  can  be  written  as  follows: 

b4" 

M  =  (1/2  +  e)  bLx  +  ebL2  -  (1/2-e)  bL3  -  -p--g-  —  .  (19) 
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Converting  to  coefficient  form  by  dividing  by  l/2p V?  ( 2b)  ( 1 )  ( 2b) 
yields 


CM  =  i/2 ( 1/2  +  e)  CL 


+  !  CL, 


-  1/2 ( 1/2-e)  C 
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APPENDIX  II 

WIND  TUNNEL  WALL  CORRECTIONS 


Wind  tunnel  boundary  corrections  were  developed  by  using  a 
computer  program  provided  by  NASA  (based  on  the  formulation  in 
References  9  and  16)  to  calculate  C^,  C^,  and  phase  for  an  air¬ 
foil  oscillating  in  a  wind  tunnel.  The  data  tables  in  Reference 
17 were  then  used  for  comparative  free-air  values.  The  theoret¬ 
ical  values  were  based  on  compressible  thin-airfoil  theory. 

The  method  of  images  was  used  to  simulate  the  walls. 

The  magnitude  ratios  of  CN  and  C^  in  free  air  to  the  values  in 
a  tunnel,  as  well  as  phase  differences,  are  shown  in  Figures  52 
and  53  for  pitching  and  plunging  motions.  The  corrections  on 
magnitude  are  largest  at  low  values  of  k  and  high  Mach  number. 
The  corrections  are  all  less  than  10  percent  except  for  the 
pitching  moment  in  translation,  where  corrections  up  to  25  per¬ 
cent  apply  at  small  k.  The  values  of  CM  for  that  condition, 
however,  are  in  the  0.002  range  and  are  therefore  close  to  the 
experimental  resolution  threshold.  This  makes  the  correction 
meaningless.  Phase  differences  between  test  and  theory  are  less 
than  10  degrees,  except  for  translation  where  the  correction  is 
15  degrees. 

The  corrections  calculated  here  were  found  to  be  small  in  com¬ 
parison  with  the  large  charges  in  both  magnitude  and  phase  due 
to  stall;  therefore,  they  have  not  been  included  in  the  wind 
tunnel  results. 
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APPENDIX  III 

STATIC  (NONOSCILLATORY)  TESTS 


Each  airfoil  was  tested  under  steady-a  condition  over  the 
entire  range  of  a's  and  M's  for  which  dynamic  data  were  obtained. 


Figures  54  through  56  show  static  Cjg  and  Cm  variations  with  a 
for  both  airfoils.  The  lines  faired  through  the  test  points  on 
the  Cn  plots  are  drawn  to  show  an  abrupt  drop  at  asTALL  for 
M  =  0.2  and  0.4  to  conform  to  the  trends  of  the  dynamic  data 
and  to  other  Boeing  experience  with  the  same  airfoil  sections. 


The  increased  values  of  CNmaX  an<^  “STALL  shown  by  the  cambered 
sectior  at  M  =  0.2  and  0.4  are  consistent  with  results  reported 
elsewhere,  such  as  Reference  13.  At  M  =  0.6,  the  Cn  curves  of 
both  airfoils  are  dominated  by  strong  compressibility  effects 
over  which  camber  appears  to  have  little  influence. 


Both  profiles  show  small  values  of  pitching  moment  about  the 
quarter-chord  at  a's  below  those  for  stall.  The  data  for  M  = 
0.2  show  considerable  irregularity  in  the  form  of  the  curves. 
This  is  due  to  the  low  dynamic  pressure  and  the  .  ulting  low 
transducer  signal  level.  The  Cm's  shown  for  the  0012  airfoil 
at  M  =  0.2  and  0.4  near  zero  lift  are  slightly  negative, 
although  zero  would  have  been  expected  because  of  symmetry. 

The  deviation  is  attributed  to  a  zero  drift  of  the  transducers. 
Note,  however,  that  the  stall  phenomena  of  principal  interest 
in  this  test  program  generally  involve  Cm's  which  are  orders 
of  magnitude  larger.  Consequently,  the  discrepancy  is  unim¬ 
portant.  Steady-flow  drag  data  for  the  23010-1.58  airfoil  are 
shown  in  Figure  55.  No  drag  data  were  measured  for  the  NACA 
0012  airfoil. 


The  static  characteristics  of  the  23010-1.58  airfoil  with  3- 
degree  trailing-edge  reflex  are  shown  in  Figure  57  for  M  =  0.4. 
There  is  no  perceptible  difference  in  the  Cn  data.  The  low  -a 
Cm  data  are  shifted  from  a  slightly  negative  value  to  a  con¬ 
sistent  zero  level.  While  the  change  is  within  the  range  of  the 
drift  mentioned  above,  it  also  agrees  with  previous  results 
such  as  those  reported  by  Davenport^. 
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PITCHING  MOMENT  COEFFICIENT  <Cm)  NORMAL  FORCE  COEFFICIENT  (CN) 


Figure  54.  Static  CN  and  CM  Characteristics  for  the 
Vertol  23010-1.58  Airfoil. 


Figure  55.  Steady  Drag  Data  for  Vertol  23010-1.58 

Airfoil  From  Wake  Traverse  Measurements. 
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APPENDIX  IV 

DISCUSSION  OF  INSTANTANEOUS  WORK 

The  instantaneous  work  function  6W  is  shown  in  Figures  58 
through  60  for  the  Vertol  23010-1.58  airfoil  as  a  function  of  o. 
Positive  6W  indicates  work  done  on  the  airfoil  by  the  tunnel 
flow  and  can  lead  to  instability  if  the  system  is  free  to  re¬ 
spond.  It  is  important  to  note  that  the  area  inside  the  6W-a 
curve  is  not  the  integrated  work  function  AW  (or  damping)  shown 
on  the  output  listings.  The  SW  function  must  be  integrated  with 
respect  to  time  rather  than  to  a  to  obtain  AW.  The  incremental 
work  function  does,  however,  show  the  amount  of  work  input  or 
output  at  each  angle  of  attack  during  a  cycle  of  oscillation. 
Since  all  data  have  been  averaged  over  five  cycles  or  more,  the 
instantaneous  work  function  presented  here  represents  an 
averaged  value. 

Figure  58  shows  the  effect  of  increases  in  a0  on  the  6W  func¬ 
tion.  Since  the  7.3-  and  9.6-degree  a0  data  do  not  stall 
during  the  cycle,  the  loop  is  small,  with  a  somewhat  larger 
area  below  the  ordinate  indicating  stability.  At  a„  =  12.2 
degrees,  the  airfoil  stalls  at  15  degrees,  causing  a  nose-down 
pitching  moment.  Since  a  is  still  increasing  up  to  17.5  degrees, 
6W  is  negative.  When  da/dt  changes  sign  at  17.5  degrees,  6W 
bt  ’omes  positive  since  Cm  and  da/dt  are  now  both  negative.  This 
is  the  condition  when  aerodynamic  forces  can  drive  the  airfoil. 

At  a„  of  14.9  and  17.5  degrees,  the  negative  6W  is  increased 
in  a  duration  and  in  magnitude,  since  the  6W  break  still  occurs 
at  15  degrees.  This  shows  that,  for  a  constant  frequency,  the 
6W  and  therefore  the  Cm  break  (since  da/dt  does  not  change 
rapidly  at  15  degrees)  are  independent  of  ae  to  a  first-order 
approximation.  The  final  loop  on  Figure  58  at  a„  =  24.6  degrees 
shows  work  input  on  the  negative  da/dt  portion  and  work  output 
during  the  positive  portion  of  the  cycle.  The  net  integrated 
value  of  work  per  cycle  (AW  or  damping)  is  shown  to  be  stable 
for  this  case  in  Figure  17. 

Figure  59  shows  the  effect  of  frequency  on  6W.  Ths  positive 
increment  due  to  stall  is  seen  to  be  largest  when  the  stall 
occurs  at  the  top  of  the  a  excursion.  This  occurs  at  49  Hertz 
for  this  a0  (12.5  degrees). 

The  effect  of  a  decrease  in  Mach  number  from  0.4  to  0.2  is 
shown  in  Figure  60  for  two  frequencies,  16  and  96  Hertz,  for 
a0  near  stall. 

The  sudden  increase  in  6W  seen  earlier  does  not  occur  until 
approximately  22.5  degrees  for  low  trequencies,  as  seen  from 
the  first  two  loops  in  Figure  60.  At  the  next  higher  a0  of  24.6 
degrees,  the  airfoil  is  fully  stalled  over  the  whole  cycle.  For 


the  high-frequency  oscillation,  5W  is  mostly  negative  over  the 
whole  cycle  at  a0  of  15  and  17.5  degrees.  At  24.6  degrees, 
both  negative  and  positive  values  are  present.  The  net  integra¬ 
ted  AW  is  unstable  as  shown  in\Figure  16.  Figure  61  shows  the 
effect  of  Mach  0.6  for  16  and  96  Hertz.  For  the  low-frequency 
oscillation,  the  6W  values  are  negative  below  an  a0  of  7.5 
degrees  for  both  positive  and  negative  pitch  rates. 

Above  7.5  degrees,  the  loop  is  essentially  symmetrical  about 
the  ordinate  axis  and  therefore  indicates  that  the  cycle  damp¬ 
ing  is  neutral  or  positive.  For  high-frequency  oscillation, 

6W  is  negative  over  three-quarters  of  the  cycle,  with  the 
only  positive  inputs  occurring  between  maximum  and  mean  a  on 
the  pitch-down. 

Figure  62  shows  6W  for  the  symmetrical  airfoil.  a0  is  increased 
from  9.9  to  14.7  degrees  in  the  top  sequence.  Stall  occurs 
2.5  degrees  earlier  for  the  symmetrical  airfoil  than  for  the 
cambered  airfoil  in  Figure  54.  The  three  figures  on  the  bottom 
show  6W  characteristics  for  a  sequence  of  increasing  frequency 
at  constant  an .  The  data  at  97  Hertz  are  more  unstable  than 
the  corresponding  cambered  airfoil  data  in  Figure  59. 

From  these  incremental  work  plots,  it  can  be  seen  that  the 
damping  per  cycle  can  be  a  very  misleading  parameter,  since  the 
instantaneous  work  or  force  inputs  can  be  much  larger  than  the 
cycle  average.  The  damping  per  cycle  can,  however,  be  used  to 
predict  the  limit  oscillation  when  such  an  oscillation  is 
established. 
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Figure  61.  Variation  of  Incremental  Work  for  the 
Vertol  23010-1.58  Airfoil  at  M  =  0.6, 
f  =  16  and  96  Hertz,  and  Aa  =  5°. 
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Incremental  Work  for  the 
1.58  Airfoil  at  M  =  0.6, 
Hertz,  and  Aa  =  5°. 
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Figure  62.  Variation  of  Incremental  Work  for  NACA  0012 
(Modified)  Airfoil  at  M  =  0.4. 
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An  experimental  investigation  into  rotor  blade  dynamic  stall  was 
conducted.  Two  typical  airfoils,  the  NACA  0012  (modified)  and  the 
Vertol  23010-1.58,  were  tested  by  measuring  the  differential 
pressures  acting  during  pitching  and  translatory  oscillations  over 
a  range  of  angles  of  attack  covering  the  stall  regime.  There  were 
two  main  objectives.  The  first  was  to  determine  the  influence  of 
stall  on  the  aerodynamic  pitching  moment  of  an  airfoil  executing 
pitching  motions  corresponding  to  the  elastic  torsional  oscillations 
of  a  rotor  blade.  The  second  was  to  determine  the  extent  to  which 
time- varying  angle  of  attack  could  affect  the  maximum  aerodynamic 
lift  that  a  rotor  blade  can  develop.  The  cambered  profile  remained 
free  of  negative  damping  for  3  to  6  degrees  beyond  the  point  where 
the  symmetrical  section  encountered  negative  damping  in  the  Mach 
0.2  to  0.4  range.  The  cambered  section  exhibited  higher  attainable 
CN's  than  the  symmetrical  section  under  all  conditions  of  motion, 
except  for  the  highest  frequ*.  cy  tested  at  Mach  0.6. 
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